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EVALUATION 
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Sperry  Research  Center  was  successful  in  obtaining  the 
smooth  impulse  response  for  aircraft  like  targets.  This 
result  has  application  for  computing  the  backscatter  or 
bistatic  response  for  any  type  aircraft  at  any  aspect  angle. 
The  advantage  of  this  solution  is  that  the  response  is  over 
the  whole  radar  spectrum.  One  can  readily  compute  what  the 
response  is  for  any  radar  transmitted  waveform.  Therefore, 
these  results  have  application  over  all  types  of  radars  and 
can  be  used  for  system  evaluation.  Evaluation  of  applying 
radar  absorbing  material  has  shown  that  8-10  db  reductions 
can  be  accomplished  over  certain  frequencies.  This  technique 
allows  the  design  engineer  the  opportunity  to  choose  the  best 
frequencies  for  detection  or  reduction.  The  success  of  this 
program  has  added  to  the  Air  Force  goal  for  computing  the 
response  of  complex  targets  at  all  frequencies  to  use  in 
radar  system  development,  as  outlined  in  RADC  TPO  2E. 


DANIEL  L.  TAURONEY 
Project  Engineer 
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INTRODUCTION 

This  document  is  submitted  as  the  final  report  in  response  to 
the  requirements  set  forth  in  Contract  No.  F30602-77-C-0166  between  the 
Sperry  Corporate  Research  Center  (SCRC) , Sudbury,  Massachusetts,  and  the 
Air  Force  Systems  Command,  Rome  Air  Development  Center  (RADC) , Griffiss 
Air  Force  Base,  New  York.  The  general  objective  of  this  program  is  to 
develop  techniques  for  predicting  the  radar  response  of  airborne  targets. 

Predicting  the  wideband  response  of  targets  for  a particular  radar 
involves  two  fundamental  seeps.  The  first  is  to  calculate,  from  the  given 
structual  geometry  of  the  target,  the  target  impulse  response  or,  equiva- 
lently, the  target  frequency  response  at  all  frequencies.  This  impulse 
response  is  the  fundamental  electromagnetic  characterization  of  the  target. 

It  is  the  knowledge  of  this  impulse  response  that  enables  the  response 
for  any  incident  waveform  to  be  simply  calculated.  The  second  step  in 
predicting  the  radar  response  is  to  apply  standard  linear  systems  theory 
to  develop  a wideband,  coherent,  radar  range  equation  which  includes  the 
effects  of  sky  and  receiver  noise. 

In  this  effort  the  techniques  for  calculating  the  smoothed 
impulse  response  from  triangular  fins,  quadrilateral  fins  and  from  more 
realistic  aircraft  models  have  been  developed  and  are  described  in  Section 
2.  The  responses  are  computed  by  the  direct  time  domain  solution  of  the 
space-time  integral  equation.  Results  obtained  for  fins  and  a MIG-21  model 
using  this  technique  were  verified  with  direct  time  domain  measurements 
and  are  discussed  in  Section  2. 

At  high  frequencies  when  the  characteristic  target  dimensions 
are  many  wavelengths,  the  target  scattering  model  can  be  reduced  to  a 
set  of  returns  from  scattering  centers  on  the  target.  These  scattering 

-1- 

» 

I 


center  returns  are  due  to  regions  of  th  target  where  the  surface  deriva- 
tives are  discontinuous  and  yield  a singularity  in  the  impulse  response 
of  the  target.  Since  these  returns  result  from  local  regions  in  space-time, 
time  domain  techniques  provide  the  natural  way  to  solve  for  them.  In  Sec- 
tion 3 time  domain  techniques  are  developed  to  calculate  the  impulse 
response  of  these  scattering  centers.  Again  the  space-time  integral  equa- 
tion approach  provides  a very  efficient  and  accurate  technique  for  these 
calculations.  Results  are  presented  for  the  prolate  spheroid,  right 
circular  cylinder  and  flat  plate  triangles  as  a function  of  aspect  angle. 

For  scattering  centers  separated  by  several  wavelengths  the  total  response 
is  then  a superposition  (in  the  time  domain)  of  individual  responses. 

The  time  domain  measurement  facility  is  described  in  Section  4 
and  smoothed  impulse  response  results  are  displayed  for  triangular  fins 
and  MIG-21  models.  Measurements  of  the  scattering  center  responses  of  a 
large  triangular  fin  and  a large  circular  cylinder  are  also  presented. 

These  scattering  center  responses  were  measured  for  both  conducting  targets 
and  radar  absorbing  material  coated  target. 

A computer  method  is  developed  in  Section  5 for  calculating  the 
broadband  radar  response  from  the  target  impulse  response.  This  method 
is  both  highly  efficient  and  accurate  and  provides  the  noise-computed, 
pulse-compressed  video  response  of  these  targets  when  illuminated  by  a 
wideband  radar.  A model  of  the  ALCOR  radar  is  used  to  demonstrate  these 
techniques . 


SECTION  2 


SPACE-TIME  INTEGRAL  EQUATION  SOLUTION  FOR  REFINED  AIRCRAFT  MODEL 

2.1  GENERAL  PROBLEM 

The  space-time  integral  equation  (STIE)  is  valid  for  any  shape 
of  conducting  target.  Solution  techniques  have  been  developed  for  a 
number  of  important  classes,  such  as  capped  and  flat-end  cylinders,  pro- 
late spheroids,  wires,  flat  plates,  and  others  [1-4].  In  Ref.  1,  the 
response  of  a simple  aircraft  model  with  rectangular  wings  was  computed. 
In  the  present  effort,  the  aim  was  to  treat  a more  realistic  aircraft 
model:  the  aircraft  chosen  for  analysis  was  a MIG-21,  illustrated  in 

Fig.  1.  The  model  for  calculation  purposes,  illustrated  in  Fig.  2,  con- 
sists of  quadrilateral  wings  (basically  triangular)  on  a thin  cylindrical 
fuselage.  Central  to  this  problem  is  the  solution  of  the  STIE  on  trian- 
gular and  general  quadrilateral  plates.  The  solution  to  this  problem 
and  the  extension  to  plates  on  a cylinder  is  presented  in  this  section. 
Comparison  with  measurements  are  also  made  in  this  section. 

All  targets  considered  here  have  plane  symmetry;  the  direction 
of  incidence  is  in  the  plane  of  symmetry;  and  the  polarization  is  taken 
to  be  perpendicular  to  this  plane  (TE  mode).  (See  Fig.  2.)  These 
restrictions  resuit  in  more  efficient  computation,  and  permit  the  compu- 
tational results  to  be  directly  verified  by  time  domain  measurements. 

It  should  be  noted  that  these  restrictions  are  not  inherent  to  the  STIE 
technique  and  can  be  removed  at  the  expense  of  greater  computational 
complexity . 

The  incident  field  is  taken  as  a smoothed  impulse,  which  yields 
a response  valid  from  DC  to  well  into  the  resonance  region.  The  tech- 
niques of  Section  3 (direct  computation  of  scattering  center  responses) 
coupled  with  the  high-frequency-augmentation  technique  described  in 
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Ref.  5 can  then  be  used  to  obtain  the  impulse  response,  or  equivalently, 
the  frequency  response  over  the  entire  spectrum. 


2.2  STIE  FORMULATION 

In  the  STIE  approach  to  the  scattering  problem,  the  solution 
is  found  directly  in  the  time  domain  by  a procedure  of  stepping  in  time. 
The  surface  currents  induced  on  the  target  by  the  incident  field  and  the 
mutual  interactions  of  the  surface  elements  are  computed  first  for  all 
time,  using  a set  of  integro-dif ferential  equations  (the  STIE),  subject 
to  the  appropriate  boundary  conditions.  Once  the  surface  currents  have 
been  found,  the  scattered  field  in  any  direction  can  be  calculated,  using 
another  STIE,  by  considering  the  equivalent  problem  of  Fig.  3.  In  this 
equivalent  statement  of  the  problem,  the  conducting  surface  has  been 
replaced  by  source  currents  radiating  into  free  space. 

The  equations  describing  the  surface  solution,  contain  terms 
which  reflect: 

(a)  the  direct  influence  of  the  incident  field, 

(b)  the  influence  of  cylinder  currents  on  other  cylinder 
currents , 

(c)  the  influence  of  fin  currents  on  other  fin  currents, 

(d)  the  influence  of  cylinder  currents  on  fin  currents, 

(e)  the  influence  of  fin  currents  on  cylinder  currents. 

The  derivation  and  solution  of  these  equations  has  been  des- 
cribed in  detail  in  the  above  references  and  will  therefore  be  presented 
here  only  in  outline.  The  solution  to  the  general  quadrilateral  plate, 
not  treated  elsewhere,  will  be  discussed  in  greater  detail. 

2.2.1  Incident  Field 

Ideally,  one  would  choose  an  impulse  as  the  time  domain  incident 
field,  since  a simple  convolution  of  the  impulse  response  with  an  actual 
incident  waveform  would  yield  the  actual  response.  Moreover,  the  impulse 
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FIG.  3 Equivalent  problem  of  scattering  from  surfaces  with  fins  attached. 
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response  is  intimately  related  to  the  actual  geometry  of  the  target, 
suggesting  the  possibility  for  techniques  that  determine  the  impulse 
response  by  inspection  of  geometry. 

For  the  numerical  solution,  however,  it  is  not  practical  to 
use  an  ideal  impulse  for  excitation.  Instead  a regularized  or  smoothed 
impulse  is  used.  The  form  of  this  illumination  at  the  origin  is  the 
Gaussian  regularization  of  an  impulse,  namely 


H1(t) 


n 


2 2 
-n  t 


(1) 


which  has  unit  integral  over  all  time  and  converges  to  a delta  functional 
as  n goes  to  infinity.  The  time  domain  integral  equations  can  be  solved 
for  bodies  with  linear  dimensions  up  to  several  pulsewidths  of  this 
smoothed  impulse.  In  this  solution  consideration  is  limited  to  bodies 
of  this  size. 

The  vector  H is  assumed  to  lie  in  the  plane  of  symmetry 
(x  = 0).  The  direction  of  incidence  makes  an  angle  a with  the  z-axis. 

(See  Fig.  2.) 

Let  r be  the  unit  vector  pointing  to  the  source,  and  $ the 

-►i 

angle  of  E with  respect  to  the  x-axis.  Then 


e E1 


1 

c 


f xH1 


For  time  in  units  of  light  meters  and  for  the  geometry  of  Fig.  2,  in 
cartesian  coordinates  at  the  point  r: 


P = 


r • 


A 

r 

o 


r = (0  , sin  a , cos  a) 
o 
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( 2 2 

H(r,t)  * e n ^t  + P^  (sin  8,  - cos  a cos  8,  sin  a cos  8) 

/ /fr 


9e 

e 3t  (r,t) 


-n  (t  + p)  | -2-  e n (t  + | (cos  8,  cos  a sin  8,  -sin  a sin  8) . 


2.2.2  Farfield 


An  integral  equation  for  the  scattered  field  may  be  obtained 
by  considering  the  equivalent  problem  of  Fig.  3,  in  which  the  scatterer 
has  been  replaced  by  source  currents  J which  are  equal  to 


J = a x H , 
n 


by  virtue  of  the  conducting  surface  boundary  condition.  Using  the  free 
sPace  Green  s function , the  solution  to  the  EM  wave  equation  then  becomes 


H(r,t)  = 


■i‘'“  •*/!&•**]“•• 


T)  x aR  } dS* 


where 


T = the  retar  ;ime  t - R/c 
-► 

r * position  vector  to  the  observation  point 

r'  = position  /ector  to  the  integration  point 

■ -►  -+  . 

R = |r  - r'| 


R R 

c peed  of  light 


a = normal  to  the  surface. 

n 

In  most  applications,  the  quantity  of  interest  is  the  far  scat- 
tered field.  In  the  limit  that  )r|  becomes  large  the  first  term  in  the 
integral  becomes  negligible. 


a 

r 


1 

R 


1 

r 


We  can  thus  write 


•+s 
r H 


(r,t) 


-L  f 

471  J 


9 J (r* , 1 ) 
3t 


x a dS' 


t - R 


(4) 


with  time  expressed  in  light  meters. 

2.2.3  Surface  Equations 

To  obtain  the  integral  equation  for  currents  on  the  surface  of 
a conducting  body  with  finite  thickness,  we  can  start  from  equation  (3) . 
The  boundary  conditions  are  that 


a x H 
n 


and  that  the  electromagnetic  field  is  zero  inside  the  body.  Moving  the 
observation  point  r to  the  surface  of  the  conductor,  applying  a limiting 
procedure , we  obtain 


J(r,t) 


2an  x H (r,t) 


with  time  expressed  in  light  meters. 
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This  equation  is  used  to  solve  for  the  surface  current  on  the 
thick  part  of  the  target,  i.e.  the  cylinder.  Since  J under  the  integral 
is  evaluated  at  the  retarded  time,  T,  this  equation  can  be  solved  numeri- 
cally by  stepping  in  time.  Note  that  the  integral  is  over  the  entire 
surface  and  thus  contains  contributions  from  the  fins  as  well  as  from 
the  rest  of  the  cylinder.  This  formulation  does  not  apply  to  open  thin 
surfaces,  since  the  boundary  condition  that  H = 0 on  the  other  side  of 
the  surface  is  not  satisfied. 

A thin  surface  is  one  in  which  E and  H do  not  vary  across  its 
thickness.  In  our  equivalent  problem  of  Fig.  3,  the  thin  surface  is  thus 
replaced  by  a current  sheet.  The  boundary  condition  for  a conducting 
surface  is  that  the  tangential  component  of  the  total  E is  zero,  hence 


E = [ 
tan 


■ LES  . S‘] 


tan 


on  the  surface 


(6) 


The  total  electric  field  in  space  can  be  written  in  terms  of  the  magnetic 

-► 

and  electric  potentials  (A  and  <J>)  : 


E(r,t)  = Ei(r,t)  - p " A ^ Lt'  “ V <J>  (r,t)  , 


where,  if  A and  <f>  are  related  by  the  Lorentz  gauge: 


(7) 


y . ^ + e ii  =0 

dt 


(8) 


then 


A(r 


^jdS’ 

T = t - R 


(9) 


Equations  (7)  and  (8)  combined  with  the  surface  boundary  con- 
dition (6)  yields 
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(10) 


[»('  • «.«)  - - -[' *^1 

Jtan  tan 

where  tan  denotes  the  tangential  component. 

Replacing  A by  equation  (9)  results  in  a first  kind  integro- 
differential  equation,  where  the  unknown  J appears  only  under  the  integral 
sign.  Again,  note  that  the  (retarded)  integral  is  over  the  entire  surface 
and  represents  the  influence  of  cylinder  as  well  as  fin  currents. 


Equation  (10)  is  to  be  solved  subject  to  the  fin  edge  conditions. 


~r-  -r 

These  conditions  are  known  for  J at  edge  point  r^: 


J, (r  , t)  = 0 
1 e 


9J1  - 

“ST  re " ° 


J ( 
II 


r.t]  - /? 


(free-space  edge) 


(11) 


J (r  ,t)  = 0 (fin-cylinder  join) , 
II  e 


where 


d is  the  perpendicular  distance  from  the  edge 

1 and  II  refer  £o  the  perpendicular  and  parallel 
components  of  J. 

The  join  edge  conditions  also  hold  at  the  symmetry  plane,  for  the  case 
of  a fin  by  itself. 

Unlike  the  second  kind  integral  equation  for  the  cylinder  body 
(equation  5) , equation  (10)  requires  the  evaluation  of  second  order  space 
derivatives.  It  is  this  aspect  that  makes  this  equation  more  difficult 
to  solve  numerically. 
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2.3  NUMERICAL  SOLUTION  OF  ST IE 


2.3.1  General 

The  space-time  integral  equations  (5)  for  the  cylinder  and  (9) 
and  (10)  for  the  fins  are  solved  numerically  by  representing  the  surfaces 
as  approximately  equal  patches.  The  current  is  assumed  to  be  constant 
over  a patch  at  a particular  time,  except  for  edge  patches  where  the  anal- 
ytic variation  is  known.  The  centers  of  patches  are  called  "nodes".  The 
distance  factor  R is  also  assumed  to  be  constant  over  the  integration 

patch,  except  for  the  "self-patch"  which  contains  the  observation  node 
-► 
r . 


Consider  first  equation  (9)  on  the  fin.  We  separate  "self" 
and  "non-self"  contributions  to  the  integral,  obtaining: 


A(ri,t)  = Yi>J(ri,t)  + Ans  (r^  ,t) 


(12) 


where 


J(r . ,T) 


'WV*1  ’ « E-S2—  “j 


(13) 


The  proportionality  factor  Y^  is  obtained  by  analytic  integration  over 
the  self-patch  and  is  given  by 


Yi 


AS. 

l 


a purely  geometrical  quantity. 


(14) 


Hence  we  can  write  the  surface  current  as 


'(vw  ■ t;  (‘‘v11  - v,;i’el) 


(15) 


-r- 

The  term  Aj^g  is  found  by  the  retarded  integration  of  equation  (13)  , while 
is  found  by  solution  of  the  differential  equation  (10).  This  solution 


f 
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is  pursued  further  in  the  next  subsection. 


! 


To  write  equation  (5)  for  a cylinder  node  numerically,  we  proceed 
similarly  by  separating  out  the  self-term,  obtaining 


“ 2Sni  X H1(Vfc)  +|Pi  5{'i't)]|  + 


(16) 


where 


-sEv  «sJ“j 


(17) 


The  meaning  of  the  term 


11 


is  as  follows:  Let  a , a be  the  direc- 

u v 

tions  of  principal  curvature  on  the  surface  with  curvatures  K ,K  . Then 

u v 


|[Vl]|  = 'i(S»3„  'V1  - av3v  <V>) 


(18) 


where 


/AS  1 /K  - K \ 

'i  ■ J /-r  (JLrJ)l 

This  is  derived  in  detail  in  an  appendix  to  Ref.  1.  We  cam  thus  solve 
for  J in  terms  of  its  components 

VV"  • (ih)  jwV>  * 2(\iSi<v‘>  S»)J 

VV>  - (it?)  {wv‘>  * Oj  • 

The  current  at  a cylinder  node  is  thus  given  directly  by  the  retarded 
integration  over  the  surface  and  the  incident  field. 


(19) 


(20) 
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Coordinate  transformations  are  necessary  to  evaluate  (20) : 

■+ 

from  the  Cartesian  system  in  which  the  J are  stored  to  the  local  u,v  sys- 

tern,  and  back  again.  (Alternatively,  J could  be  stored  in  the  u,v  system, 

but  then  coordinate  trans formations  are  required  to  evaluate  I „ . ) 

NS 

. 

J and  A are  computed  at  discrete  time  points,  separated  by  At. 

To  evaluate  J(T)  and  its  first  time  derivative,  an  interpolation  is  per- 
formed on  the  five  points  in  time  nearest  T.  These  points  are  smoothed 
linearly  and  a quadratic  is  fit  through  the  resulting  three  points.  The 
process  of  smoothing  reduces  the  problem  of  numerical  instabilities  in 
the  solution. 

2.3.2  Detailed  Solution  of  STXE  on  Fin 

We  now  discuss  in  greater  detail  the  solution  of  equations  (9) , 
(10)  on  the  flat  plate,  subject  to  the  edge  conditions  (11).  In  our  geom- 
etry, the  flat  fin  is  located  in  the  plane  y = 0 (see  Fig.  2).  Then  ex- 
panding equation  (10)  into  components 


92A  92A  32A  32A 

x _ x z y 

2 3x3z  3x3y  + C 3t 


3t 


32a 


3t 


3x 


32a 


32a 


32a 


3z 


3E1 

x 


3e 


5.  + * + 1 + E — 5. 

2 3z3x  3z3y  3t 


(21) 


Using  a finite  difference  approach,  these  equations  are  solved  by  forming 
the  numerical  space  derivatives  of  Afr^.t)  based  on  the  values  of  A at 
surrounding  nodes.  Representing  the  time  derivative  also  in  the  finite 
difference  formulation 


32A 

3t2 


(t) 


A (t + At)  - 2A(t)  + A(t  - At) 
(At)2 


(22) 
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we  can  thus  solve  for  A(t  + At)  , based  on  the  values  of  A at  (t)  and 
(t-At)  , on  the  numerical  space  derivatives  of  A at  (t)  , and  on  the  inci- 


dent field:  e.g.  for  the  x-component 


A (t  + At)  = 2A  (t)  - A (t-At)  + (At)' 

X XX 


92A  (t)  92a  (t)  92a  (t)  9E1(t) 

X z y X 

a 2 9x9z  9x9y  9t 


all  evaluated  at  the  nodes  r^  on  the  fin.  Two  problems  present  themselves: 

1.  How  to  form  the  space  derivatives,  so  that  a stable 
solution  results; 

2.  The  numerical  differentiation  is  invalid  at  and  near 
the  edges,  where  the  edge  conditions  have  to  be  in- 
voked to  derive  a value  for  A. 

“►  -► 

The  edge  conditions,  however,  are  known  for  J and  not  for  A.  On  the  other 
hand,  A is  bounded  at  the  edge,  as  will  be  shown,  so  that  A may  be  obtained 
directly  by  integration  of  equation  (9) , making  use  of  the  known  variation 
of  J. 

Again  we  separate  the  self-patch  from  the  integration: 

*(Vt}  = + Vvt}  ' (24 

where  A^g,  given  by  equation  (13)  contains  the  contribution  from  all 
patches  (including  the  cylinder  and  other  fins) , except  the  self-edge 
patch  at  re>  and  Ag  is  the  integral  over  the  self  patch: 


vV>  - i J lsiiP£l 


The  details  of  this  last  integration  are  left  for  Appendix  8.1.  The 
result  is  that 
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ASII(Vt) 


(25) 


where  II  and  1 refer  to  the  parallel  and  perpendicular  components.  The 
proportionality  constants  y , y^  are  purely  geometrical,  but  reflect  the 
variation  of  J (equation  11)  over  the  patch. 


Y„  = 0 


1 2 AS  ' [2AS 

3 e 

8^  IT  ' 'l  16%  TT 

I2AS ~ 

’ ~ TT 

-r 

'l  “ 48%  tt 


at  a free  edge; 


at  a bound  or 
symmetry  edge; 


at  a free  comer; 


Yu  =0 


1 4AS 

y±  = 2 — i 

8^  TT 


at  a bound  corner. 


(The  directions  II , 1 of  the  bound  corner  are  defined  to  be  those  of  the 

bound  edge  attached  to  it.)  The  free  corner  constants  are  derived  from 

the  assumption  that  J + 0 at  the  corner,  and  that  the  corner  is  about 

90°.  A more  obtuse  corner  could  be  treated  as  an  edge.  The  bound  comer 

constants  are  obtained  from  the  observation  that  one  edge  is  bound  and 

the  other  free.  (Confusion  is  avoided,  when  it  is  noted  that  for  all 

these  constants  y as  well  as  for  (19)  and  (14)  , the  quantity  under  / 

represents  p or  the  radius  of  a circular  patch  with  the  area  As  . ) If 

e 

M is  the  transformation  from  Cartesian  to  the  edge  (#,i)  coordinate 
system. 


J 

II 

= M 

z 

J, 

J 

L XJ 

X 
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then 


+ + r -I  +ef  f 

As(re't)  = 0]  J <re'fc) 


(27) 


where 


[y]=M 


-1 


Y»  ° 


0 Y, 


This  value  of  A can  then  be  used  in  (24)  to  obtain  A at  the  edge. 

S 

-►p  ff  , 

The  effective  current,  J , is  also  found  by  integration  over 
the  edge  patch  in  accordance  with  the  edge  conditions  (see  Appendix  8.1). 
The  result  is 


J«ff=k  J 
II  I II  0 


JCff  =*  k J , 
i 1 10 


where 


2.226  — ■ — i , k = .4244 

(2AS  /IT)* 
e 


/ 


2As  /IT 


at  a free  edge; 


= 0 , k^  = 1 at  a bound  or  symmetry  edge; 


(4AS  /it)  * 

k = k = 1.2113  at  a free  corner; 

11  1 ^d 


k = 0 , k = 2.226 

II  l 


A 


(4AS  A) 
e 


at  a bound  corner. 


(29) 


d is  the  perpendicular  distance  from  the  edge  node  to  the 


point  r ; 
o 
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J = J(r  ) , obtained  by  interpolating  between  the  two  nearest  interior 
o o 

nodes : 

Jq  = ai  Jfr^)  + a2  J(^2)  . (30) 

The  edge  geometry  is  illustrated  in  Fig.  4.  Again  we  need  to  apply  the 
coordinate  transformation  M,  obtaining 

-►e  f f -*■  r « i -►  -► 

J (re,t)  = [6]  JQ(ro,t)  , (31) 

where 

[6]  - H-1 

We  have  thus  obtained  the  effective  current  at  the  edges  based  on  interior 
currents  and  the  edge  conditions.  These  effective  values  are  used  in 
forming  the  integrals.  The  values  of  A at  the  edge  are  used  for  evalua- 
ting  space  derivatives  of  A at  interior  nodes. 

Another  detail  is  required  for  the  evaluation  of  equation  (23) : 

Since  observation  nodes  exist  only  on  the  fin  plane  which  has  y = 0,  it 

2 

is  not  possible  to  evaluate  9 A^/9x9y  numerically.  Instead  we  perform 
9/9x  numerically  and  obtain  9A^/9y  from  the  integral  equation. 


_9_ 

9y 


AyU.t) 


J_J_  f Jy  (r ' ,T) 

4TT  9y  J R 


dS' 


_9_ 

9y 


A 

y 


( r , t) 


ds' 


(32) 


Note  that  nodes  on  the  fin  itself  do  not  contribute  to  this  term  on 
account  of  the  factor  (R  • y) , which  is  zero  when  observation  and  source 
point  are  both  in  the  y = 0 plane. 
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2.3.3  Summary  of  Solution  Procedure  on  Fin 


We  can  now  summarize  the  solution  procedure.  At  time  t: 


A.  First  compute  at  interior  nodes  r^;  assuming  that  A(r^,t) 


has  been  computed  by  procedure  at  previous  time  step 


1.  Obtain  A^ (r\  ,t)  by  integration  of  J(r',T)  over  all 


other  patches,  via  equation  (13) 

*►  ”►  ■> 

2.  Solve  for  J(r^,t)  using  A,  A^s  and  equation  (15). 


B .  Compute  at  edge  nodes  r 


1.  Obtain  A (r  ,t)  by  integration  of  J(r',T)  over  all 
NS  e 


other  patches,  via  equation  (13). 


2.  Find  J from  results  of  step  (A2)  and  equation  (30). 
o 


-fg  ff  “► 

3.  Solve  for  J (r  ,t)  using  J and  equation  (31) . 

e o 


-►  -►e  f f 

Solve  for  A,(r  ,t)  using  J and  equation  (27). 
o e 


Solve  for  A(rg,t)  using  Ag , A^s  and  equation  (24). 


We  now  have  A and  J at  all  nodes  of  the  fin  up  to  time  t. 


C.  Solve  differential  equation  at  interior  nodes  r^ . 


1.  Evaluate  9A^/9y  (r,t)  by  equation  (32). 


Evaluate  space  derivatives  of  A(r^,t)  using  values  of 


A at  surrounding  nodes. 


3.  Evaluate  incident  fields  at  (r^,t) 


4. 


Obtain  A(r^,  t + At)  from  equation  (23). 


We  are  now  ready  to  repeat  the  process  at  the  next  time  step. 
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Notes : 


1.  This  procedure  is  performed  after  calculation  of  surface 
currents  on  cylinder  at  time  t. 

2.  Interaction  of  fins  and  cylinder  occurs  via  integrals  (17) 
and  (13)  which  are  over  all  patches. 

3.  Computational  storage  and  time  is  reduced  by  invoking 
symmetry:  surface  currents  are  not  stored  or  computed 
at  image  points  (in  the  present  case,  points  x < 0) , but 
their  contributions  are  included  in  the  integrals. 

2.3.4  Surface  Derivatives  on  Flat  Surface 

It  remains  to  discuss  the  problem  of  forming  the  derivatives 
numerically  on  the  surface  of  the  flat  fins.  In  earlier  formulations  of 
the  thin  plate  problem  [3] , a rectangular  grid  was  used  (Fig.  5a) . The 
derivatives  were  then  evaluated  by  standard  finite  difference  formulas. 

For  example, 

2-*-  ■*  ■*  -*■ 
d A j+j  = A(r  +■  Ax)  +■  Air -Ax)  - 2A(r)  . 

3x2  (Ax) 2 

This  method  requires  modification  when  the  plate  outline  is  not  rectangular 
and  the  grid  cannot  be  made  to  coincide  with  the  plate  edges. 

A skew  hexagonal  grid  is  illustrated  in  Fig.  5b.  This  grid  is 
still  regular  in  the  sense  that  all  intervals  Ax  are  equal  and  all  intervals 
Az  are  equal.  The  angle  0^  is  restricted  to  <45°  (for  larger  values  of 
0^,  the  grid  can  be  rotated  by  90°). 

A skew  coordinate  system  x' , z'  is  defined  (see  Fig.  5b)  and  the 
derivates  are  found  by  the  chain  rule.  For  example 
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FIG.  5 Grid  choices. 


f 


I 


where 


32A 

(r) 

A(r  + Ax' ) 

+ A(r  - Ax' ) 

- 2A(r) 

3x'2 

(Ax')2 

32A 

<r) 

A(r  + Az' ) 

+ A(r  - Az' ) 

- 2A(r) 

3z'2 

(Az’)2 

and  where 

i ■+  Ax*  \ 1 /*^  ■+  Ax*  \ 

A(r  + Az' ) = i Ap  + Az'  + j A Ir  + Az'  + ^-> 

with 

Az'  = Az  f Ax'  = Ax/cos  0^ 

Although  still  not  completely  general  (the  outside  dimensions 
have  to  be  commensurate  with  the  Ax  and  Az  which  in  turn  are  related  to 
each  other  by  0^  and  02) , this  grid  permits  a large  class  of  targets  to 
be  treated.  Most  aircraft  wings,  for  example,  can  be  represented  closely 
by  this  configuration.  Variations  of  the  illustrated  case  are  the  trap- 
ezoid, parallelogram,  and  triangle.  The  hexagonal  grid  concept  was 
tested  by  comparing  the  calculation  of  a square  with  hexagonal  grid  with 
the  results  for  a rectangular  grid.  The  latter  had  previously  been  com- 
pared with  measurement  (3).  Although  some  differences  were  observed  in 
the  intermediate  calculation  of  surface  currents,  the  scattered  farfield 
results  were  indistinguishable.  This  response  is  illustrated  in  Fig.  6; 
it  is  for  a 2x2  plate  with  incident  pulse  n = 1 at  normal  incidence. 

The  skew  hexagonal  grid  was  used  for  all  the  response  calcula- 
tions presented  in  the  next  section. 

While  the  regular  hexagonal  grid  was  adequate  for  present 
purposes,  an  effort  was  made  to  generalize  the  representation  of  surface 
derivatives  in  order  to  make  the  solution  procedure  truly  independent  of 
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plate  geometry.  One  method  was  to  write  a second  order  Taylor  expansion 
*►  *► 

of  A about  the  point  r.  The  coefficients  (i.e.,  the  derivatives  at  r)  were 

then  found  by  a least  square  fit  of  this  polynomial  to  the  values  of  A at 

surrounding  nodes.  Such  a procedure  would  allow  an  arbitrary  arrangement 

of  nodes.  Although  the  derivatives  obtained  this  way  had  reasonable 

values,  the  insertion  of  this  technique  in  our  time-stepping  STIE  solution 

caused  numerical  instabilities  which  increased  in  time.  Moreover,  these 

instabilities  were  worst  the  more  asymmetrically  the  surrounding  nodes  were 

placed.  Even  for  a completely  regular  grid,  however,  the  solution  was 

unstable.  The  reason  for  this  result  is  apparently  as  follows:  the  fitting 

procedure  assumes  constant  values  of  second  order  derivatives  over  the 

entire  region  of  the  fit.  (In  contrast,  the  simple  finite  difference 

method  assumes  only  the  constancy  of  one  second  derivative  along  each  of 

2 2 

the  axes  of  evaluation;  e.g.  formula  33  assumes  only  that  8 A/8x  is  con- 
stant along  the  x-axis.)  For  our  problem  this  is  a bad  approximation, 
which  effectively  results  in  discontinuities  in  second  and  first  derivatives 
in  going  from  one  patch  to  an  adjacent  patch. 

A solution  to  this  problem  may  be  a higher  order  development 
of  A so  that  derivatives  can  be  made  continuous  across  patch  boundaries. 

This  is  precisely  what  is  done  in  the  procedure  for  solving  differential 
equations  known  as  the  finite  element  technique.  This  technique  was 
investigated.  It  was  decided  that  the  complexity  of  implementation  did 
not  warrant  its  use  in  the  present  problem.  However,  for  more  general 
problems,  it  appears  to  be  the  method  of  choice  and  it  is  therefore 
presented  in  outline  in  Appendix  8.2. 

2.4  COMPUTATIONAL  RESULTS 

We  have  already  presented  in  Fig.  6 the  results  of  plate  response 
calculations  for  a square  plate  with  a hexagonal  grid,  which  showed  good 
agreement  with  previous  calculations  made  with  a normal  finite  difference 
grid.  As  a further  test  of  the  method,  calculations  were  made  for  45° 
and  30°  right  triangles  (standing  on  am  image  plane)  as  a function  of 
incidence  angle.  Finally,  the  response  of  two  triamgles  attached  to  a 
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right  circular  cylinder,  representing  the  model  of  a MIG-21,  was  computed. 
All  these  calculations  were  verified  by  measurement  on  the  time-domain 
scattering  range. 

2.4.1  Triangle  Responses 

The  45°  triangle  has  4”  sides.  The  30°  triangle  has  sides  of 
4"  x 6.9282"  with  the  long  side  on  the  symmetry  plane.  Note  that  the 
free-space  equivalent  of  this  geometry  is  the  triangle  augmented  by  its 
image . 


The  response  of  the  45°  triangle  is  given  in  Figures  7 and  8 
for  a at  45°  intervals  from  0°  to  180°.  The  polarization  is  TE.  The 
incident  smoothed  impulse  has  n = 0.5 (/inch);  that  is,  the  pulsewidth, 
given  between  2%  points  as  W = 4/n,  is  8".  The  pulsewidth  used  in  the 
measurements  was  slightly  shorter  than  this.  The  measured  response  was 
transformed  by  a convolution  procedure  to  the  same  (n)  so  that  direct 
comparison  could  be  made.  The  transformed  measured  resp'  es  are  super- 
posed in  dotted  lines  on  the  calculated  response  curves. 

The  response  of  the  30°  triangle  is  given  in  Figures  9 and  10. 
Other  parameters  are  the  same  as  for  the  45°  triangle.  These  responses 
are  discussed  further  in  Section  4. 

In  Figure  11  the  effect  of  shape  changes  is  illustrated.  Res- 
ponses are  given  for  variations  on  the  shape  of  the  small  fin  of  the  MIG 
model.  All  fins  have  the  same  area.  The  fin  illustrations  are  drawn  to 
the  same  scale  as  the  time  axis  (ct)  of  the  responses.  Also  shown  to 
scale  is  the  incident  pulse.  The  effect  of  the  shape  change  is  not 
pronounced.  It  is  seen  that  the  solution  for  trapezoidal  shape  responses 
have  a problem,  in  that  there  is  residual  oscillation.  For  this  reason, 
the  MIG  model  responses  were  calculated  with  triangular  wings.  The  part 
of  the  aircraft  response  due  to  the  large  fins  should  not  be  affected  by 
this  approximation  since  the  wing  is  very  nearly  triangular;  the  part  of 
the  response  due  to  the  small  fin  may  be  affected  somewhat  (according  to 
Fig.  11) , but  the  magnitude  of  the  small  fin  response  is  small  in 


comparison  with  the  cylinder  or  large  fin  responses. 


2.4.2  Fins  on  Cylinder 

Figure  12  illustrates  the  responses  at  broadside  incidence  of 
the  small  fin  by  itself,  of  the  cylinder  representing  the  fuselage,  and  of 
the  small  fin  on  the  cylinder.  The  pulsewidth  parameter  here  is  n=l. 
(Hence,  the  2%  pulsewidth  is  4.)  The  cylinder  length  is  14.8  with  a radius 
of  0.81.  The  small  fin  is  a 30°  triangle  with  height  1.73.  These  dimen- 
sions, if  taken  as  inches,  represent  a 1/32  model  of  the  MIG  fuselage  and 
stabilizers . 

The  cylinder  broadside  specular  return  looks  somewhat  like  a 
doublet  centered  at  t = -1.62,  with  its  backswing  modified  due  to  nature 
of  the  curved  specular  surface.  (See  Section  3 for  analysis  of  cylinder 
response.)  The  last,  positive,  peak  is  the  return  from  the  creep  wave, 
occurring  at  about  t = 2.5.  The  fin  by  itself  has  a response  somewhat 
like  a doublet  centered  near  t = 0,  followed  by  a creep  return  near  t=3. 

It  has  the  effect  on  the  response  of  cylinder  + fin  of  modifying  the  nega- 
tive swing  of  the  cylinder  response.  In  the  total  response,  the  creep 
return  of  the  cylinder  is  still  present  but  is  modified  and  followed  by 
another  echo  at  t = 6 , which  is  interpreted  as  cylinder-fin  interaction, 
being  due  to  a creep  from  fin  to  fin  around  the  backside  of  the  cylinder. 
Ringing  is  seen  beyond  t = 10;  however,  the  calculation  suffers  from 
numerical  noise  beyond  this  time. 

Figure  13  illustrates  the  responses  at  broadside  incidence  of 
the  large  fin  by  itself  and  on  the  cylinder.  The  fin  is  a 30°  triangle 
with  height  4.  The  response  of  this  fin  is  similar  to  that  of  the  small 
fin  except  for  larger  amplitude,  different  relative  peaks,  and  creep  return 
time  due  to  the  difference  in  relative  length  of  fin  and  incident  pulse. 

The  magnitude  of  the  fin  response  is  about  the  same  as  the  specular  return 
of  the  cylinder  and  hence  causes  a large  modification  around  t = 0 of  the 
cylinder  response  when  the  fin  is  added  to  the  cylinder.  The  creep  response 
of  the  cylinder  is  still  evident  near  t = 3,  and  is  followed  by  fin-cylinder 
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creep  return  near  t = 8.  The  total  response  of  the  MIG  model  is  shown  in 
Fig.  14,  also  at  broadside  incidence.  At  this  aspect,  to  first  order, 
this  looks  similar  to  the  large-fin-on-cylinder  response,  the  main  effect 
of  the  additional  fin  being  to  increase  the  fin  part  of  the  response.  At 
the  time  of  the  creep  return,  there  is  a noticeable  effect  due  to  fin-fin 
interaction  in  addition  to  fin-cylinder  interactions.  Beyond  t = 10  or  so, 
the  results  are  probably  invalid  due  to  numerical  noise. 

The  response  of  the  MIG  model  was  measured.  In  order  to  permit 
comparison,  both  measured  and  calculated  responses  were  processed  by  a 
convolution  procedure  to  yield  a response  for  pulsewidth  parameter  n = 0.5. 
The  result  of  this  is  shown  in  Fig.  14b.  The  agreement  is  good  for  the 
initial  positive  and  negative  swings.  There  appears  to  be  a relative  time 
shift  in  the  responses  after  t = 5,  which  is  when  the  result  of  fin-fin 
and  fin-cylinder  interaction  occurs.  The  discrepancy  is  not  understood, 
but  may  have  something  to  do  with  the  detailed  treatment  of  fin  currents 
at  the  singular  edges. 
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SECTION  3 

STIE  CALCULATION  OF  SCATTERING  CENTER  RESPONSES 


The  principal  contributors  to  the  response  of  a complex  target 
are  the  scattering  centers  (surface  discontinuities  and  points  on  a sur- 
face with  a normal  bisecting  the  angle  between  incidence  and  observation) . 
This  is  more  true  at  higher  frequency,  when  the  separation  between  scatter- 
ing centers  is  greater  than  a few  wavelengths,  or  equivalently,  when  the 
separation  is  larger  than  the  incident  puiselength,  and  the  direct  inter- 
action between  scattering  centers  via  surface  currents  (creep  waves)  can 
be  neglected.  This  fact  has  permitted  the  successful  application  to 
response  calculations  of  physical  optics  theory  (which  holds  only  in  the 
limit  of  very  high  frequency  and  fails  to  predict  polarization  dependence) 
and  the  geometrical  theory  of  diffraction,  GTD  (which  is  based  on  the 
response  of  a wedge  and  includes  polarization  dependence,  but  requires  the 
application  of  a transformation  from  wedge  to  actual  scattering  center 
geometry) . 

In  the  time  domain,  the  impulse  response  predicted  by  these 
theories  can  be  viewed  as  a series  of  singularities.  These  response  sin- 
gularities can  also  be  obtained  analytically  by  considering  the  surface 
currents  in  the  immediate  neighborhood  of  the  scattering  center.  (See 
Ref.  1 and  Appendix  8.3).  Such  an  approach  yields  the  physical  optics 
terms  and  a first  order  (polarization  dependent)  correction.  This  method 
is  made  use  of  in  what  follows  as  an  approximate  verification  of  the  new 
technique. 

The  time  domain  (STIE)  technique  presented  here  is  a direct 
method  of  calculating  the  broadband  response  of  scattering  centers,  edges 
as  well  as  points  on  smooth  surfaces.  The  response  of  a complex  target 
can  then  be  obtained  by  direct  superposition  in  the  time  domain.  Alter- 
natively , these  scattering  center  responses  can  be  used  as  input  to  the 
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"high  frequency  augmentation  technique"  [5]  together  with  the  smoothed 
impulse  response  as  calculated  in  Section  2,  to  yield  the  complete  impulse 
response . 

♦ 

3 . 1 METHOD 

The  approach  is  the  well  tested  technique  of  the  space-time 
integral  equation  described  in  Section  2,  but  using  as  the  incident  field 
a smoothed  impulse  of  very  short  duration,  T.  The  far field  response  of 
the  scattering  center  is  then  considered  to  occupy  a time  window  equal  to 
the  pulselength  and  corresponding  to  the  time  that  the  scattering  center 
is  illuminated  (Fig.  15) . The  region  of  the  target  around  the  center 
which  contributes  to  this  response,  are  then  those  points  which  are  within 
distance  cT/2  along  the  direction  of  incidence.  (We  consider  here  only 
the  monostatic  case,  although  the  method  is  readily  applicable  to  the 
bistatic  case  as  well.) 


In  the  numerical  implementation,  the  surface  and  time  are 
divided  into  a fine  mesh.  The  surface  currents  are  calculated  first  and 
then  from  them,  the  farfield.  (See  Section  2.2.)  Considerable  efficiency 
of  calculation  is  realized  when  use  is  made  of  the  fact  that  the  space-time 
region  on  the  surface  is  bounded  by  the  lines  p = ct  and  p = cT  - ct 
(see  Fig.  15d) , where  p is  the  projected  distance  of  the  surface  node  from 
the  scattering  center  along  the  direction  of  incidence. 

The  method  is  demonstrated  on  a prolate  spheroid  and  a right 

circular  cylinder  as  a function  of  aspect  angle.  Each  of  these  objects 

is  centered  at  the  origin  with  axis  along  the  z-axis.  The  direction  of 

incidence  is  in  the  yz  plane  at  an  angle,  a,  with  respect  to  the  z-axis. 

In  the  response  curves,  the  time  of  return,  T , is  measured  with  respect 

z 

to  a return  from  the  origin  and  is  thus  given  by  cT^  = -2pc  where  p^  is 
the  distance  of  the  specular  point  from  the  origin  along  the  direction  of 
incidence.  The  response  to  both  TE  and  TM  polarizations  are  presented, 
where  TE  means  E transverse  to  the  yz  plane.  The  incident  pulsewidth 
parameter  was  n = 16.  The  pulsewidth  for  these  calculations  was  T = 8/n. 
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FIG.  15  Scattering  center  response  calculation  relationships. 


The  calculations  were  verified  by  measurement  using  a large  target 
and  short  pulse  excitation. 


The  prolate  spheroid  responses  as  a function  of  a,  essentially 
give  the  result  due  to  a smooth  scattering  center  with  principal  curvatures 
K^,  K2>  The  right-circular  cylinder  responses  are  the  responses  of  a curved 
90°  edge  (K^  = 00 , = 1/a)  as  a function  of  aspect.  In  addition,  the  end-on 

(a  = 0)  and  broadside  (a  = 90°)  responses  of  the  right-circular  cylinder 
represent  the  specular  responses  of  a flat  disc  (K^  = K2  = 0)  and  a line 
(K^  = 1/a,  K9  =0)  on  a curved  surface. 

3.2  RESULTS 

3.2.1  Prolate  Spheroid 


A typical  response  of  the  prolate  spheroid  is  illustrated  in 
Fig.  16;  this  one  is  for  TM,  a = 60°.  At  other  aspects,  the  responses 
are  visually  dully  similar  so  that  it  is  more  instructive  to  present  the 
results  in  tabular  form. 


The  principal  curvatures  along  the  directions  a and  a are 


K 


1 


£1  r 
2 2 ' *2 
a b 


2 

a 


(34) 


where 

„ 2.2  2 2 

B = a sm  a + b cos  a 

a,b  are  the  semi -minor  and  -major  axes, 
a^  is  along  a meridian  line, 

a2  = $.  (See  Fig.  15e) 

The  expected  impulse  response  for  a smooth  surface  with  finite 
curvatures  is  of  the  form  (see  Appendix  8.3) 
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— h(— ) = P S(T)  + S U (T)  + R T U ( T)  + 
a V a / z z z 


normalized  to  r^/a  = observation  distance/target  radius,  where  P^,  , 

P are  the  impulse,  step,  and  ramp  coefficients,  and 
z 


(35) 


„ cT 
ct  z 


I 


T = 


All  terms,  except  the  first  are  polarization  dependent.  This  expression 

was  fitted  to  the  computed  responses  in  order  to  facilitate  interpretation. 

The  results  are  presented  in  Table  I.  Also  shown  in  the  table  are  the 

theoretical  values  for  P and  S as  obtained  from  the  first  order  correc- 

z z 

tion  to  the  physical  optics  response  (Appendix  8.3).  For  most  cases,  the 

agreement  is  seen  to  be  close.  The  coefficient  R is  not  very  accurate 

z 

due  to  the  arbitrary  nature  of  the  second  order  fit.  Also,  it  is  seen 

that  the  impulse  coefficient  P shows  a false  sliaht  polarization  depen- 

z 

dence  - this  is  due  to  numerical  error.  The  P term  is  due  to  the  zero- 

z 

order  component  of  the  physical  optics  response.  As  can  be  seen  in  Fig. 

16,  this  term  dominates  the  calculated  response.  The  accuracy  of  the 

calculation  could  probably  be  improved  by  computing  ''correction  currents" 

on  the  surface,  i.e.  the  total  current  minus  the  physical  optics  currents. 

-*  A -*inc 

' i the  surface  STIE  of  Section  2.2.3,  the  latter  is  the  term  J_„  = 2a  x H 

po  n 

.2.2  Right  Circular  Cylinder 

The  geometry  for  these  cases  is  illustrated  in  Fig.  15f.  The 
end-on  and  broadside  aspect  responses  are  presented  in  Fig.  17.  The  end-on 
response  is  polarization  independent  and  is  given  theoretically  by 


— h(— ) = D 6'  (T)  + P 6 (T)  + H I-**  u 
a \ a / z z z 


(T)  + 


(36) 
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FIG.  17  RCC  scattering  center  responses  - end-on  and  broadside. 
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where  D , P are  the  doublet  and  impulse  coefficients, 
z z 


H is  the  coefficient  of  the 
z 


functional , 


The  expected  values  (from  physical  optics)  of  the  first  term  is  D =0.50. 

z 

At  broadside  the  response  is  polarization  dependent  and  is  given  by 
(Appendix  8.3) 


-3/2  -1/2 

G T ' u(T)+H  T u(T)  + ... 

z z 


(37) 


the  expected  values  are  G = .15915,  H = -.079  (TE) , H = +.079  (TM) . 

z z z 

The  results  of  fitting  these  functionals  to  the  measured  responses  are 
presented  in  Table  II.  Agreement  with  theoretical  expectations  is  good 
except  for  the  90°  TE  case,  which  is  probably  in  error. 


At  intermediate  angles  the  responses  are  much  smaller  in  magni- 
tude than  at  0°  or  90°.  Responses  are  plotted  at  the  same  scale  in 
Fig.  18  for  both  TE  and  TM  cases. 

II 

The  theoretical  form  of  the  response  is  given  by 


— h/— I = H T-^  u ( T)  + K u(T)  + ...  (38) 

a \ a / z z 

The  results  of  fitting  these  functionals  to  the  calculated  responses  is 
given  in  Table  II. 

Measurements  of  a RCC  are  described  in  Section  4.  The  measure- 
ments were  made  with  a pulselength  about  3 times  longer  than  these  calcula- 
tions, so  that  comparison  (after  convolution)  can  only  be  made  over  the 
initial  portion  of  the  measured  response.  The  results  agree  within  this 
time  window. 
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TABLE  II 

RCC  Scattering  Center  Responses 


Aspect 

Polarization  cT  /a 

D 

G 

P 

H 

K. 

z 

z 

z 

z 

z 

z 

0 

NA  -2.0 

.4930 

0. 

. 162 

-.051 

30 

TE  -2.7321 

0. 

0. 

0. 

. 3725 

-.63 

30 

TM 

0. 

0. 

0. 

.2067 

-.004 

45 

TE  -2.8284 

0. 

0. 

0. 

.2930 

-.34 

45 

TM 

0. 

0. 

0. 

.1492 

.07 

60 

TE  -2.7321 

0. 

0. 

0. 

.2823 

-.39 

60 

TM 

0. 

0. 

0. 

.1574 

.02 

90 

TE  -2.0 

0. 

- .1585 

0. 

.006 

-.07 

90 

TM 

0. 

- .1599 

0. 

.073 

-.29 

SECTION  4 

TIME  DOMAIN  MEASUREMENTS 

A necessary  ingredient  for  any  successful  development  of  tech- 
niques for  the  solution  of  previously  unsolved  scattering  problems  is  the 
availability  of  actual  measurements  that  can  be  used  to  verify  the  new 
techniques.  During  this  effort  the  Sperry  Research  Center's  Time  Domain 
Scattering  Range  was  used  for  these  measurements.  Section  4.1  describes, 
in  summary  form,  this  measurement  system,  its  operation,  and  its  important 
performance  characteristics.  The  measurement  results  which  were  taken 
for  this  effort  are  described  in  Section  4.2. 

4.1  MEASUREMENT  SYSTEM 

The  Sperry  Research  Center's  Time  Domain  Scattering  Range  is  a 
system  for  obtaining  a low-noise  response  in  the  nanosecond  region.  It 
consists  of  a ground  plane  scattering  range,  a step  function  source,  a 
sampling  oscilloscope  receiver,  and  a laboratory  instrumentation  computer 
for  control  and  processing.  The  system  signal  source  is  a high-voltage 
switch  which  generates  a 300  V step  function  with  a risetime  less  than 
100  ps . The  signal  is  radiated,  virtually  undistorted,  from  a wire  trans- 
mitting antenna  protruding  through  a circular  ground  plane  20  feet  in 
diameter.  This  wave  illuminates  the  target  and  the  resulting  scattered 
waveform  is  received  on  a coaxial  horn  antenna,  which  essentially  smoothes 
and  differentiates  the  signal  and  thus  provides  the  smoothed  impulse 
response  of  the  target.  The  received  waveform  is  sampled  by  a 12  GHz 
oscilloscope  that  has  been  triggered  by  the  initial  pulse  and  whose  sam- 
pling gate  deflection  is  under  the  control  of  a small  instrumentation 
computer.  Unprocessed  data  are  displayed  on  the  oscilloscope  CRT  while 
the  sampled-and-held  waveform  is  passed  through  a low-pass  filter,  digitized, 
read  into  the  computer,  and  stored  on  magnetic  tape  automatically.  This 
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system  has  been  designed  to  correct  the  long-term  timing  drift  and/or 
amplifier  drift.  In  addition,  the  waveforms  are  stored  in  such  a way  that 
they  are  ready  for  the  subsequent  operations  of  averaging  (to  remove  short- 
term noise)  and  baseline  processing.  The  effects  of  a time  varying  baseline 
are  subtracted  from  measured  waveforms  to  improve  system  accuracy. 

The  salient  characteristics  of  the  range  are  the  speed  and  sim- 
plicity with  which  multi-octave  frequency-domain  data  can  be  obtained. 

These  advantages  accrue  because  the  time-domain  scattering  range  yields  an 
"uncontaminated"  interval  of  time  between  the  arrival  of  the  direct  wave 
and  the  arrival  of  unwanted  reflections.  Targets  are  usually  located 
anywhere  from  two  to  five  feet  from  the  transmitting  antenna.  The  response 
from  the  antenna  tip  and  the  table  edge  occur  at  approximately  15  ns.  Thus, 
a "clean  window"  exists  between  4 ns  and  15  ns  which  can  be  used  to  view 
the  target  responses.  The  entire  region  between  the  direct  transmission 
and  the  table  edge  response  forms  a convenient  time  "window"  to  view  the 
target  response  and  allows  one  to  "gate  out"  (in  time)  unwanted  reflections. 
Thus,  undistorted  transient  target  responses  can  be  viewed  without  resorting 
to  elaborate  and  expensive  anechoic  chambers.  In  addition,  a single  time- 
domain  measurement  obviates  the  requirement  for  tedious  measurement  of  the 
amplitude  and  phase  responses  at  many  frequencies. 

The  accuracy  of  the  measurement  system  has  been  estimated  for 
the  results  presented  in  this  section.  The  peak  of  the  incident  pulse 
as  measured  on  the  sampling  oscilloscope  is  approximately  400  mV,  and  a 
typical  target  response  has  a peak  value  in  the  vicinity  of  10  mV.  When 
using  the  10  mV  scale  on  the  sampling  oscilloscope,  the  standard  deviation 
of  the  sample  mean  is  estimated  to  be 

0 = 0.5  mV 

v 

if  16  scans  are  averaged.  Thus,  the  estimated  standard  deviation  of  the 
sample  mean  v is  in  the  vicinity  of  5%  of  the  peak  value  of  the  target 
response . 
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In  addition,  the  measured  responses  are  further  processed  by  means 


of  a convolution  procedure  to  obtain  the  response  due  to  a Gaussian  shaped 
incident  pulse  rather  than  the  approximate  smoothed  impulse  used  in  the 
actual  measurements.  Figure  19  displays  the  actual  measured  incident  pulse 
along  with  the  smoothed  Gaussian  pulse  that  was  used  in  the  convolution 
process.  The  frequency  spectrum  of  these  two  pulses  is  displayed  in  Fig. 
20.  The  time  domain  expression  for  the  Gaussian  pulse  is  given  by 

e(t)  = E — exp  f-a2  (t  - t ) 2 1 

° /F  L n ° J 

and  the  frequency  domain  expression  is 


E (to)  = Eq  exp|-  M2/(2an)2]  exp(-joJto)  . 

The  amplitude  coefficient  Eq  was  set  equal  to  the  dc  value  of  the  measured 

incident  pulse.  The  width  coefficient  a was  obtained  by  requiring  that 

n 

the  50%  value  of  |e(w) | occur  at  the  same  point  in  frequency  as  the  50% 
value  of  the  magnitude  of  the  transformed  measured  incident  pulse.  The 
width  of  the  resulting  Gaussian  pulse  becomes  0.59  ns  or  approximately  7 
inches.  An  added  benefit  of  this  process  is  the  reduction  of  high  fre- 
quency noise  in  the  response. 

This  measurement  system  provides  the  capability  of  measuring  and 
recording  the  processed  high  precision  smoothed  impulse  response  data 
on  computer  computable  tape.  These  data  may  then  be  used  to  find  the 
response  due  to  any  radar  waveform  whose  spectrum  is  contained  within  the 
spectrum  of  the  original  measurement  smoothed  impulse  illumination.  Further 
details  of  this  system  may  be  found  in  the  references  [1-7]. 

4.2  MEASURED  RESULTS 

In  previous  studies  [1-6]  the  smoothed  impulse  response  was 
measured  for  numerous  target  geometries . These  previous  measurements 
are  summarized  in  Table  3.  During  tne  present  contract,  smoothed  impulse 
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(a)  MEASURED 


(b)  GAUSSIAN 


ra-2M 


FREQUENCY  (GHz) 
(b)  GAUSSIAN 


FIG.  20  Spectrum  of  incident  pulse. 


TABLE  III 


Summary  of  Previous  Measured  Responses 


June  1970  [2]  and  May  1973  [5] 


ADC 

satellite 

model 

UES 

satellite 

model 

GGTS-2 

satellite 

model 

Wire 

L/10  = 10 

SSS 

satellite 

model 

EWRS 

satellite 

model 

November  1970  [6)  and  May  1973  [5] 


8 in.  cube 

Right  square  cylinder,  L = 12  in.,  D = 4 in. 

Right  circular  cylinder,  L = 12  in.,  D = 4 in. 

Right  circular  cylinder,  L = 8 in.,  D = 4 in. 

4i  in.  sphere 
8 in.  sphere 

Sphere  capped  cylinder,  L = 12  in.,  D = 4 in. 
Sphere  cone  sphere 


May  1974  [31 


8x8  in. 
4 x 8 in. 
8 in. 
4x8  in. 
8x8  in. 
8x8  in. 
8 in. 


square  plate 
square  plate 
circular  disc 
corner  reflector 
circular  cylinder  section 
parabolic  cylinder  section 
parabolic  dish 


March  1976  [7] 


Cylinder  with  square  fins,  L 
Cylinder  with  square  fins,  L 


8 in. , D = 4 in. 
12  in. , D = 4 in. 
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responses  were  measured  for  a number  of  triangular  fins,  a Revelle  MIG-21 
model,  and  a machined  MIG-21  model.  In  addition  smoothed  impulse  scattering 
center  responses  were  measured  using  a large  triangular  fin  and  a large 
right  circular  drum.  The  responses  were  measured  for  both  the  naked  con- 
ducting bodies  and  for  Radar  Absorbing  Material  (RAM)  covered  bodies.  The 
specific  measurements  taken  on  this  effort  are  summarized  in  Table  4. 

The  Revelle  MIG-21  model  is  1/32  scale  and  is  displayed  in  Fig.  21(a). 

The  machined  MIG-21  model  that  was  used  for  the  measurements  is  shown  in 
Fig.  21(b)  and  the  dimensions  are  given  in  Fig.  2.  Note  that  the  photos 
in  Fig.  21  were  taken  of  the  actual  models  sitting  on  a mirror  surface. 

Fig.  22  displays  the  actual  models  sitting  on  the  time  domain  scattering 
range.  The  target  geometries  of  the  triangles  is  shown  in  Fig.  23.  The 
target  geometries  of  the  large  triangle  and  the  large  drum  are  displayed 
in  Fig.  24.  The  smoothed  impulse  responses  which  are  displayed  in  this 
section  are  due  to  the  incident  smoothed  impulse  shown  in  Fig.  19(b) , which 
has  a width  of  0.59  ns. 

In  Fig.  25,  the  measured  smoothed  impulse  response  of  a MIG-21 
Revelle  model  is  shown  for  the  case  where  the  image  plane  is  perpendicular 
to  the  wings  and  coincides  with  the  rudder.  This  model  was  fabricated  from 
Revelle  kit  H267,  cut  in  half,  and  painted  with  Emerson  & Cuming,  Inc. 
Eccoshield  ES  conducting  paint.  This  provides  the  response  for  the  case 
of  TE  polarization.  For  0°  incidence  (nose-on)  the  return  from  the  nose 
which  should  appear  at  1.2  ns  is  embedded  in  the  system  noise.  At  approxi- 
mately 3.1  ns  a negative  pulse  appears  that  is  due  to  the  scattering  from 
the  rear  edge  of  the  wings.  Note  that  the  return  from  the  nose,  the 
fuselage  sides  and  the  leading  edge  of  the  wings  is  negligible  relative 
to  the  return  from  the  trailing  edge  of  the  wing.  Finally,  at  approximately 
3.9  ns  a second  negative  return  appears  that  can  be  attributed  to  the 
trailing  edge  of  the  stabilizers.  As  the  angle  of  incidence  increases  from 
0°  to  90°  the  response  shape  changes  smoothly  reflecting  the  change  in 
relative  distances  of  the  scattering  centers.  At  90°  the  beginning  of 
the  return  from  the  side  of  the  fuselage  appears  at  approximately  1.7  ns 
followed  immediately  by  the  large  return  from  the  wings  and  stabilizers 
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TABLE  IV 


Summary  of  Measured  Responses 


Revelle  MIG- 21  Model 

0°, 

30°  , 

60°, 

90°, 

120° 

, 150°, 

180° 

Machined  MIG-21  Model 

0°, 

30°, 

60°, 

90°, 

120° 

, 150°, 

180° 

30°  Triangle 

0°, 

30°  , 

60°, 

90°, 

120° 

, 150°, 

180° 

45°  Triangle 

0°, 

30°, 

60°  , 

90°, 

120° 

, 150°, 

180° 

60°  Triangle 

0°, 

30°, 

60°  , 

90°, 

120° 

, 150°, 

180° 

Large  30°  Triangle 

0°, 

30°  , 

60°  , 

90°  , 

120° 

, 135°, 

150°  , 

Large  30°  RAM  Triangle 

120‘ 

135°,  150° , 

165°, 

180° 

Large  Cylinder 

0°, 

30°  , 

45°, 

60°, 

75°, 

90° 

Large  RAM  Cylinder 

30°  , 

, 45° 

, 60° 

, 75° 

, 90° 
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(a)  REVELLE  MODEL 


(b)  MACHINED  MODEL 


/8  ?98 


FIG.  21  MIG-21  models  on  a mirror  plane. 
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FIG.  23  Target  geometry  of  triangular  plates. 
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FIG.  24  Target  geometry  of  large  triangular  plate  and  large  cylinder  models. 
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that  approximate  a smoothed  doublet.  This  is  followed  by  an  oscillating 
waveform  due  to  currents  traveling  along  the  surface  of  the  wings.  Con- 
tinuing on  to  180°,  the  return  from  the  end  of  the  fuselage  and  the  rear 
of  the  stabilizers  appear  first  at  approximately  0.8  ns.  This  is  followed 
at  approximately  1.6  ns  by  a second  positive  pulse  due  to  the  return  from 
the  trailing  edge  of  the  wings.  The  response  oscillates  several  times 
before  damping  out  due  to  currents  flowing  on  the  wings.  The  return  from 
the  far  end  (the  nose  for  this  aspect)  which  should  appear  at  approx. 

3.3  ns  is  buried  in  the  return  from  these  traveling  wing  currents. 

Fig.  26  displays  the  measured  smoothed  impulse  response  of  the 
machined  MIG-21  model  (see  Fig.  2)  that  represents  the  actual  geometry 
used  for  the  space-time  integral  equation  computations  that  were  described 
in  Section  2.  The  responses  measured  for  the  machined  model  in  Fig.  26 
agree  very  closely  with  those  in  Fig.  25  measured  for  the  Revelle  model. 

The  wave  shape,  amplitude,  and  relative  spacings  of  features  are  in  good 
agreement.  The  machined  model  exhibits  a noticeable  return  from  its  nose 
for  aspect  angles  of  0°  to  90°,  whereas,  the  nose  return  for  the  Revelle 
model  is  smaller.  From  these  measurements,  the  location  of  the  major 
scattering  centers  may  be  identified  to  be  the  trailing  edge  of  the  wings, 
the  trailing  edge  of  the  stabilizers,  the  side  of  the  fuselage  and  the 
ends  of  the  fuselage. 

In  Fig.  27,  the  measured  smoothed  impulse  response  of  a triangular 
plate  with  a 30°  half  angle  as  shown  in  Fig.  23.  The  image  plane  is 
normal  to  the  face  of  the  plate  and  0°  incidence  corresponds  to  the  apex-on 
direction.  In  Fig.  27  one  can  see  that  at  0°  incidence  there  is  very 
little  return  from  either  the  apex  or  the  sides.  The  significant  return 
occurring  at  2.4  nsec  is  a negative  pulse  from  the  trailing  edge  of  this 
plate.  As  the  aspect  angle  increases  from  0°  to  90°  the  return  from  the 
other  portions  of  the  plate  emerge  and  at  90°  the  return  from  the  plate 
approximates  the  derivative  of  the  incident  pulse.  The  second 
positive  pulse  at  90°  aspect  is  due  to  currents  traveling  across  the  face 
of  the  plate  before  returning  energy  back  to  the  observer.  These  effects 
become  "stretched  out"  or  dispersed  as  the  aspect  angle  increases  to  180° 
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because  now  all  parts  of  the  plate  are  illuminated  at  different  times. 

In  Fig.  28,  the  measured  smoothed  impulse  response  of  a triangular 
plate  with  a 45°  half  angle  as  shown  in  Fig.  23.  The  image  plane  is  normal 
to  the  face  of  the  plate  and  0°  incidence  corresponds  to  the  apex  on 
direction.  For  this  target  the  return  from  the  front  edge  is  apparent 
even  at  0°  incidence,  in  the  form  of  a positive  pulse.  This  is  followed 
by  a negative  pulse  approximately  0.7  ns  later  from  the  rear  edge  of  the 
plate.  At  90°  incidence  the  return  approximates  a derivative  of  the  inci- 
dence pulse.  The  amplitude  in  this  case  is  smaller  than  for  the  30°  plate 
and  the  damping  more  rapid  as  expected.  Proceeding  from  90°  to  180°  yields 
a set  of  responses  very  similar  to  those  of  the  30°  plate,  however,  the 
amplitude  at  180°  is  greater  than  the  30°  plate  and  the  oscillation  is 
more  rapid  than  the  30°  plate.  This  is  expected  since  the  trailing  edge 
of  the  plate  is  moving  toward  vertical  giving  a more  abrupt  discontinuity 
to  the  surface  currents  at  an  earlier  time. 

In  Fig.  29,  the  measured  smoothed  impulse  response  of  a triangular 
plate  with  a 60°  half  angle  as  shown  in  Fig.  23.  The  image  plane  is  normal 
to  the  face  of  the  plate  and  0°  incidence  corresponds  to  the  apex-on 
direction.  For  this  target  the  return  from  the  front  edge  appears  at  1.2 
ns  for  0°  incidence,  the  return  from  the  far  edge,  a negative  pulse, 
appears  at  approximately  2.0  ns,  and  the  return  due  to  traveling  currents 
on  the  surface  appears  at  approximately  3.0  ns.  This  waveform  gradually 
changes  as  the  aspect  angle  increases.  At  90°  the  initial  return  is 
doublet  in  nature  although  the  impact  of  the  large  apex  angle  is  apparent 
in  the  asymmetry  of  the  initial  portion  of  this  waveform.  This  is  followed 
in  time  by  a damped  oscillation  due  to  vertical  traveling  currents.  As 
the  aspect  angle  increases  from  90°  to  180°  this  general  waveform  shape 
changes  only  slightly,  reflecting  the  compactness  of  this  particular 
target  geometry. 

Since  the  major  returns  from  the  MIG-21  model  were  due  to 
scattering  from  the  rear  edge  of  the  wings,  the  side  of  the  fuselage,  and 
the  end  of  the  fuselarge,  a large  model  of  a 30°  half  angle  triangular 
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plate  and  a large  model  of  a circular  cylinder  were  fabricated  as  shown 
in  Fig.  24.  The  results  of  smoothed  impulse  returns  from  the  scattering 
centers  on  these  targets  is  displayed  in  Fig.  30(a)  & (b) . The  return 
from  the  apex  of  the  30°  half  angle  triangle  is  negligible  as  shown  in 
Fig.  30(a).  As  the  aspect  angle  increases,  a return  from  the  apex  becomes 
noticeable  at  an  aspect  of  60°.  Note  that  the  return  from  the  far  edge 
of  this  plate  is  outside  the  5 ns  time  window  of  these  measurements  until 
the  aspect  angle  increases  to  60°  at  which  point  the  far  edge  return  can  be 
observed  at  a time  of  approximately  4.5  ns.  At  90°  the  large  doublet-like 
return  from  the  broadside  of  the  triangle  is  apparent.  In  Fig.  30(b)  the 
measured  returns  from  the  large  cylinder  are  displayed  where  0°  corresponds 
to  axial  incidence.  Note  that  at  0°  the  return  is  a very  symmetric  doublet- 
like waveform  that  provided  a good  approximation  to  the  derivative  of  the 
incident  pulse.  At  aspect  angle  of  30°  and  60°  the  return  at  approximately 
1.3  ns  is  the  scattering  center  return  from  the  near  edge  of  the  cylinder 

end.  At  broadside  incidence  (90°)  the  response  approximates  a smoothed 
-3/2 

t pseudo- function  as  expected. 

In  Fig.  31  the  measured  smoothed  impulse  response  of  the  large 
30°  triangle  is  displayed  for  aspect  angles  that  range  from  120°  to  180° 
(edge  on  incidence) . The  results  are  shown  for  the  metal  plate  by  itself 
and  for  the  metal  plate  when  it  was  coated  by  radar  absorbing  material 
(RAM)  from  the  vertical  edge  to  one  foot  away.  In  Fig.  31(a)  the  positive 
pulse  from  the  plate  edge  at  approximately  1.1  ns  is  followed  by  the  effect 
of  the  return  from  the  slanted  trailing  edge  at  approximately  0.3  ns 
later.  In  Fig.  31(b)  the  response  from  the  same  target  coated  with 
Emerson  & Cumings  Eccosorb  SF-2.5  RAM  is  shown.  Note  that  the  effect  of 
the  RAM  coating  in  the  time  domain  on  the  smoothed  impulse  response  of 
the  scattering  centers  is  to  reduce  the  amplitude,  to  broaden  the  response 
and  to  smooth  out  the  variation. 

As  an  example  of  the  frequency  domain  scattering  center  response, 
the  Radar  Cross  Section  (RCS)  of  the  scattering  center  for  this  target 
is  shown  in  Fig.  32  for  aspect  angles  of  150°  and  180°.  The  RCS  of  the 


FIG.  30  Smoothed  impulse  response  of  scattering  centers 
(hor.  scale:  0.5  ns/div.;  vert,  scale:  5 mV/div.). 


0 


(b)  180° 


FIG.  32  Radar  Cross  Section  (RCS)  of  large  30°  triangle  with  and 
without  RAM. 


target  without  RAM  (the  solid  curve)  is  compared  with  the  RCS  of  the  target 
with  RAM  (the  dotted  curve) . Note  that  for  150°  incidence  there  is  approxi- 
mately 8 dB  reduction  in  RCS  between  2.0  and  2.9  GHz  with  narrow  nulls 
which  are  down  approximately  16  dB  at  2.0  and  2.6  GHz.  For  180°  incidence 
there  is  approximately  10  dB  reduction  in  RCS  between  2.0  and  3.0  GHz, 
with  narrow  nulls  which  are  down  approximately  20  dB  at  2.0,  2.4,  and  3.0 
GHz. 

In  Fig.  33  the  measured  smoothed  impulse  response  of  the  large 
cylinder  is  displayed  for  aspect  angles  that  range  from  30°  to  90°  (broad- 
side incidence) . Results  are  shown  for  the  naked  cylinder  and  for  a 
cylinder  coated  with  Emerson  & Cumings  Eccosorb  SF-2.5  RAM.  The  response 
in  Fig.  33(a)  for  the  clean  target  shows  a scattering  center  return  from 
the  near  edge  at  aspect  angles  of  30,  45,  60,  and  75°.  The  positive  pulse 

returns  are  all  similar  in  size  and  shape.  At  90°,  however,  the  return 
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approximates  a smoothed  t pseudo- function , as  expected.  When  RAM  is 
added  to  the  vicinity  of  the  scattering  center,  the  response  shown  in 
Fig.  33(b)  results.  The  major  effects  are  a spreading,  smoothing,  and 
amplitude  reduction  of  the  smoothed  impulse  response  of  the  uncoated  RAM 
target. 

Fig.  34  displays  a comparison  of  the  RCS  in  the  frequency  domain 
without  and  with  RAM  coating.  For  a 60°  incidence  angle  the  RCS  is  reduced 
approximately  10  dB  between  1.2  GHz  and  3.0  GHz  by  the  RAM  coating.  Narrow 
nulls  of  approximately  20  dB  reduction  are  also  noted  at  approximately  1.5 
GHz,  2.4  GHz,  and  2.6  GHz.  In  Fig.  34(b)  the  comparison  for  broadside 
incidence  is  displayed.  A reduction  of  approximately  20  dB  is  observed 
between  2.2  GHz  and  2.9  GHz. 

It  is  obvious  from  these  measurements  that  direct  time  domain 
measurement  of  scattering  center  response  provides  a powerful  tool  in  the 
broadside  analysis  of  radar  systems. 


WITHOUT  RAM 

• • • • WITH  RAM 


VO  2.0  3.0  4.0 

FREQUENCY  (GHz) 

(a)  60° 


FIG.  34  Radar  Cross  Section  of  large  cylinder  with  and  without  RAM. 
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SECTION  5 


CALCULATION  OF  TARGET  RESPONSE  AS  SEEN  BY  WIDEBAND  RADARS 

The  impulse  response  is  the  fundamental  electromagnetic  charac- 
terization of  the  target.  Knowledge  of  the  impulse  response  enables  the 
response  for  any  incident  waveform  to  be  simply  calculated.  The  purpose 
of  this  section  is  to  describe  and  demonstrate  a technique  for  using  the 
target  impulse  responses  to  predict  the  received  video  waveforms  for  wide- 
band radars.  A simplified  (but  inclusive)  model  of  the  radar  is  developed 
and  implemented  on  a computer.  Its  use  in  obtaining  the  actual  radar 
video  response,  for  a target  whose  impulse  response  is  known,  is  demon- 
strated for  the  case  of  the  ALCOR  radar.  The  techniques  developed  in  this 
section  are  also  applicable  to  predicting  the  performance  of  coherent 
imaging  systems. 

5.1  BACKSCATTE RED  WAVEFORM 

The  inherent  linearity  associated  with  the  scattering  of  elec- 
tromagnetic fields  from  conducting  bodies  implies  that  the  target  behaves 
as  a linear  system  in  the  sense  that  the  scattered  field  can  be  represented 
as  the  convolution  of  the  incident  waveform  and  the  target  impulse  response 
as  follows: 

+°° 

r (t)  = f hT(t-T)  e (T)  dT  (39) 

J _oo 

where  r(t)  is  the  scattered  target  response,  e(t)  is  the  incident  pulse, 
and  hT<t)  is  the  target  impulse  response.  An  equivalent,  frequency-domain 
statement  of  the  relationship  between  the  incident  and  scattered  waveforms 
is  [ 81 

R(w)  = Ht  (w)  E(u»  (40) 
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where  R,  H , and  E are  the  Fourier  transforms  of  r,  hT  and  e,  respectively. 
This  relationship  clearly  shows  how  the  target  acts  as  a linear  filter 
on  the  incident  signal. 


5.2  COHERENT  RADAR  SIGNAL  PROCESSING 

Except  for  some  scaling  constants  and  delays  which  will  be  dis- 
cussed later,  equations  (39)  or  (40)  describe  the  waveform  at  the  receiving 

* 

antenna  terminals  of  the  radar.  Thus  e(t)  is  the  input  to  the  signal 
processing  section  of  the  radar. "*  The  signal  processing  sections  of  most 
modem  high-resolution  radars  are  coherent,  that  is,  they  manipulate  the 
phase  as  well  as  the  amplitude  of  the  received  waveforms.  For  example, 
a very  general  model  of  the  received  waveform  is 


r(t)  = b(t)  cosQjj^t  + 0(t)]  (41) 

where  b(t)  and  0(t)  are  the  amplitude  and  phase  modulations  of  the  carrier. 
A simple  trigonometric  identity  shows  that  (41)  is  equivalent  to 


r(t)  = £b(t)  cos  0(t)]cos 


U)Qt  - [b(t) 


sin 


6(t)]  sin  u)  t 


(42) 


which  can  be  written 


r(t)  = Re  | £b(t)  cos  0(t)  + j b(t)  sin  0(t)J 


,JV  I 


(43) 


where  Re  denotes  "real  part  of".  Equations  (42)  or  (43)  show  that  the 
quantities  in  the  rectangular  brackets  completely  describe  the  received 
waveform  in  the  sense  that,  if  they  are  known,  the  RF  waveform  for  a known 
carrier  can  be  reconstructed  by  (42). 


This  statement  is  true  assuming  the  frequency  response  of  the  receiving 
antenna  is  flat  over  the  bandwidth  of  the  scattered  waveform.  If  this 
is  not  true , then  the  frequency  response  of  the  antenna  becomes  an  addi- 
tional multiplicative  term  in  (40) . 

^ The  effects  of  an  IF  section  (if  any)  are  unimportant  to  this  discussion. 
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(44) 

I 

q ( t)  = b(t)  sin  0(t) 

are  the  so-called  in-phase  and  in-quadrature  components  of  the  waveform 
relative  to  the  carrier  [9] . 

The  p and  q components  of  the  received  waveform  can  be  obtained 
by  the  method  of  Figure  35.  The  filters  pass  the  slowly  varying  p and 
q components  but  heavily  attenuate  all  oscillations  at  much  higher  frequen- 
cies. The  two  outputs  in  Figure  35  are  often  called  the  baseband  signal, 
bipolar  video,  the  complex  modulation,  or  the  complex  envelope.  The  latter 
terminology  is  descriptive  because  (43)  can  be  written 


(t)  = Re  |[p(t)  + j q(t)] 
= Rejr(t)  e ^ 


jV! 


(45) 


where 


r(t)  = p(t)  + j q(t)  = b ( t)  e 


j 0 ( t ) 


(46) 


The  elaboration  of  the  preceeding  paragraphs  is  intended  to  show  that 

all  the  important,  or  information-bearing,  characteristics  of  typical  radar 

signals  are  contained  in  the  p and  q components  or  the  complex  envelope. 

It  is  then  to  be  expected  that  the  signal  processor  must  deal  only  with 
the  complex  envelope. 

5.3  COMPLEX  ENVELOPE  OF  THE  RETURN  SIGNAL 

A general  expression  for  the  pulse  applied  to  the  terminals  of 
the  transmitting  antenna  is 
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HMMhM 


s(t)  = a(t)  cos^u^t  + iHt)] 


where  a(t)  and  \|>(t)  are  the  amplitude  and  phase  modulations  impressed  on 
the  RF  carrier.  For  example,  for  a linear  FM  or  chirp  pulse. 


*<t)  - \ Bt2 


The  incident  waveform  in  the  vicinity  of  the  target,  e(t) , will 
be  a scaled  and  delayed  version  of  s(t).  Thus,  ignoring  for  the  moment 
some  constants  and  time  shifts,  the  return  signal  at  the  receiving  antenna 
terminals  is  given  by  (39)  or  (40) . The  p and  q components  of  this  signal 
are  simply 


p(t)  = £2  r(t)  cos  w t] 


o JLP 


q(t)  = [-2  r ( t)  sin  u>  t] 


where  the  subscript  LP  means  "low  pass  filtered  version  of"  the  quantities 
in  the  brackets  [10] . Using  some  simple  properties  of  Fourier  transforms 
and  some  careful  algebra,  it  can  be  shown  that  the  Fourier  transform  of  the 
complex  envelope  (or  p and  q components)  of  the  received  signal  is  given 


R(U»  = P(w)  + j Q(W)  = F[r(t)]  = H (u>  + <a  ) S(ld) 


where  S(w)  is  the  Fourier  transform  of  the  complex  envelope,  s(t) , of  the 
signal  and  Ht  is  the  Fourier  transform  of  the  target  impulse  response. 

Note  that  the  tilde  notation  indicates  a complex  or  baseband  function  of 
time  or  frequency.  Equation  (50)  shows  that  the  bipolar  video  return 
signal  (complex  envelope)  is  influenced  by  the  target  frequency  response 
in  a band  centered  about  the  radar  carrier  as  expected.  The  width  of  this 
band  is  the  bandwidth  of  the  transmitted  pulse.  The  returns  from  wider 
band  transmitted  pulses  will  contain  more  information  about  the  target 


-78- 


than  those  from  narrower  band  pulses. 


5.4  MATCHED  FILTERING  AND  PULSE  COMPRESSION 

The  full  range  resolution  and  detection  sensitivity  of  a radar 
are  obtained  by  appropriately  filtering  the  received  waveform.  The  matched 
filter  maximizes  the  output  signal  to  noise  ratio  (detection)  and  time 
compresses  phase  modulated  signals  (range  resolution).  Since  the  time 
(or  range)  sidelobes  in  the  matched  filter  output  are  often  unacceptably 
high,  exact  matched  filtering  is  often  not  done  in  practical  receivers. 

In  this  section,  however,  we  consider  matched  filtering  and  pulse  compres- 
sion to  be  identical  operations.  In  addition,  the  term  matched  filtering 
here  means  matched  to  the  transmitted  signal  in  white  noise,  that  is,  the 
target  is  assumed  to  be  an  ideal  sphere  whose  radius  is  much  larger  than 
the  radar  wavelength  so  that  it  is  illuminated  only  in  the  optical  region. 
Note  that  for  target  classification,  filtering  which  is  matched  to  the 
expected  return  from  each  different  class  of  extended  target  would  be 
useful . 


In  view  of  the  dominant  role  played  by  the  complex  envelope, 
it  is  perhaps  not  surprising  that  matched  filtering  can  be  done  in  the 
video  section.  Figure  36  is  a block  diagram  of  the  pulse  compression 
receiver  assumed  in  this  report.  The  frequency  response  of  the  matched 
filter  is 


H (<o)  = Ks‘(nj)  (51) 

Mr 

where  S is  the  Fourier  transform  of  the  complex  envelope  of  the  trans- 
mitted pulse  [10,  p.  245)  and  K is  a constant.  The  Fourier  transform  of 

the  compressed  complex  video,  r (t) , is  then 

c 


R (00) 
c 


KH  (u 
T 


00  ) 

o 


S(oo) 


(52) 


which  follows  from  (50)  and  (51) . 


The  last  equation  is  the  main  result  of  this  section.  The  com- 
pressed video  return  r (t)  is  obtained  by  taking  the  inverse  Fourier  trans- 

c 

form  of  (52) . Linear  video  A-scope  traces  are  then  simply 

v ( t)  = |rc(t) | (53) 

as  shown  in  Figure  36. 

The  following  subsections  describe  the  method  used  to  include 
noise  and  scale  the  A-scope  plots  so  that  they  are  physically  descriptive 
of  an  ALCOR  situation.  Several  examples  are  discussed  in  the  final  sub- 
sections. 

5.5  SCALING  VIA  THE  RADAR  EQUATION 

In  this  subsection,  the  constants  and  delays  neglected  above 
are  included  in  order  to  present  the  video  output  in  physically  descriptive 
units.  The  system  to  be  analyzed  is  shown  in  Figure  37.  This  analysis 
is  basically  a derivation  of  the  radar  range  equation  (11)  and  the  somewhat 
tedious  details  will  not  be  given. 

It  can  be  shown  that  the  complex  video  before  pulse  compression 
is  approximately  given  by 
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where  O is  the  target  RCS.  In  this  case,  the  compressed  complex  video  is 


r (t)  = 
c 


+°° 

/ 


hMF(t-T)  r(T)dT 


Using  (51)  and  (55), 


r (t) 
c 


+°o 

/ 


— oo 


S*(T  - t) 


dT 


from  which  an  application  of  the  Schwartz  inequality  (8]  gives  the  peak 
matched  filter  output 


max 

t 


(56) 


where  ET  is  the  pulse  energy  measured  at  the  transmitting  antenna  terminals. 
This  equation  allows  the  matched  filter  output  to  be  scaled  so  that  its 
peak  value  is  the  square- root  of  the  RCS  in  meters. 

5.6  INCLUSION  OF  NOISE  EFFECTS 

The  combination  of  both  the  antenna  and  the  receiver  noises  is 

assumed  to  be  measured  at  the  input  to  the  mixers  (Fig.  37)  and  is  assumed 

to  be  white  with  power  spectral  density  N /2  watts  per  Hz.  It  can  then 

o 

be  shown  (10,  p.  582]  that  the  covariance  function  of  the  noise  at  the  low 
pass  filter  outputs  is  given  by 

N ft 

R (T)  ” “ sinc(ftT)  (57) 

n tt 

where  ft  is  the  radian  cutoff  frequency  of  the  ideal  low  pass  filters. 
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From  this,  it  follows  that  the  covariance  function  of  the  complex  noise 
before  matched  filtering,  n(t) , is 

2N  ft 

R~(T)  = — “ sinc(^T)  (58) 

n 

wv ich  has  a rectangular  power  spectrum  of  level  2N  extending  from  -ft  to 

o 

+ft  radians  per  second. 

If  the  matched  filter  bandwidth  is  smaller  than  the  low  pass 
filter  bandwidth,  the  power  spectral  density  of  the  compressed  noise  is 
approx ima  te 1 y 


S„  (W)  = 2No|HMF(00)  |2 
n 

c 

= 2NqK2  | S (UJ)  | 2 (59) 

where  the  last  equality  follows  from  (51).  The  variance  of  the  compressed 
complex  noise  is  simply  the  total  area  under  the  spectrum  (59)  which  is 

Varfn  (t)l  = n (t)  n (t)  = 2N  K E_  (60) 

L c J c c o T 

where  the  overbar  denotes  statistical  average  and  E^,  is  the  pulse  energy 
at  the  transmitter  terminals 

The  peak  signal  to  noise  ratio  at  the  output  of  the  matched 
filter  for  an  ideal  target  is  then  [ll,  p.  144] 


As  with  signals,  resorting  to  complex  noise  is  merely  a convenient  way 
of  simultaneously  treating  the  quadrature  components. 
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peak  SNR  = 


t |rc(t)|: 
I n (t)  | 2 

1 c 1 


v K/ 


which  follows  from  (56)  and  (60).  This  is,  of  course,  the  noise- limited 
radar  range  equation. 

The  specific  examples  in  the  following  sections  were  computed 
using  the  inverse  Fourier  transform  of  (52).  The  amplitude  scaling  of 
the  A-scope  traces  was  done  using  (56)  so  that  the  peak  linear  video  levels 
have  the  units  of  meters  (square  root  of  RCS) . The  noise  levels  were  set 
using  equation  (61)  for  a one  square  meter  target.  The  results  in  the 
following  sections  give  this  SNR  which  is  equivalent  to  specifying  the 
range.  For  example,  a 12  dB  SNR  on  a one  square  meter  target  corresponds 
to  the  ALCOR  range  of  about  2000  km  (12] . 


5.7  EXAMPLE:  IDEALIZED  TWO-ELEMENT  TARGET 


The  following  simple  example  will  illustrate  the  ideas  of  the 
previous  subsections.  Assume  that  the  transmitted  pulse  has  the  wideband 
ALCOR  [12]  parameters  of  10  ysec  duration  and  500  MHz  bandwidth.  The  target 
is  assumed  to  be  two  spheres  separated  by  10  meters.  The  sphere  nearest 
the  radar  has  an  RCS,  o , of  9 square  meters.  The  other  scattering  center 
has  02  = 1 . 

In  order  to  properly  evaluate  the  inverse  Fourier  transform  of 
a scaled  version  of  (52)  , the  target  frequency  response  function  must  be 
obtained.  For  this  simple  target,  the  result  is 


v-v-rtS1.1*- 

T O Ro  4TT 


0)L 

- 3 T 


where  H relates  the  incident  field  to  the  scattered  field  at  reference 
T 

range  R^  and  L is  the  distance  between  the  scattering  centers.  Note  that 
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the  reference  range  cancels  out  of  the  properly  scaled  equation  (54) . 


Figure  38(a)  shows  the,  pulse  compressed,  scaled,  linear  video 
|r  (t)|  for  SNR  = 20  dB  (the  target  is  roughly  1300  km  from  ALCOR) . The 
time  axis  is  converted  to  range  by 

AR  - | (At) 

so  that  the  scattering  center  separations  are  readily  apparent.  Note  that 
at  this  range,  the  display  of  Figure  38(a)  gives  a useful  picture  of  the 
target  reflecting  structure.  Figure  38(b)  shows  the  same  target  at  10  dB 
SNR  (2100  km)  but  the  small  scattering  center  is  no  longer  distinctly 
visible  in  the  noise. 

5.8  EXAMPLE:  RIGHT  CIRCULAR  CYLINDER 

Since  ALCOR  is  primarily  used  for  space  surveillance  and  re-entry 
measurements,  a more  realistic  example  is  the  wideband  response  of  a large 
right  circular  cylinder  as  seen  by  ALCOR.  The  impulse  response  of  this 
target  is  given  in  a previous  report  [7]  and,  for  the  purposes  of  this 
example,  we  assume  that  the  target  is  a large  cylindrical  rocket  body  which 
is  20  meters  long  and  4 meters  in  diameter. 

Figure  39  shows  the  tranverse  magnetic  (TM)  and  transverse 

i 

electric  (TE)  wideband  video  response  of  this  target  at  a very  high  signal 
to  noise  ratio.  These  three  dimensional  plots  show  local  range  (-20  to 
+20  meters)  or  scattering  center  separation  horizontally  across  the  page 
with  aspect  angle  varying  in  the  vertical  direction.  The  lower  A-scope 
trace  in  each  figure  is  for  near  broadside  incidence  while  the  upper  traces 
represent  aspects  approaching  axial  incidence.  The  three  dominant  peaks 
in  each  trace  show  the  relative  locations  of  the  three  scattering  centers 
as  well  as  the  overall  target  length  projected  along  the  radar  line  of 
sight.  The  smaller  ripples  or  peaks  are  the  range  or  time  sidelobes 
characteristic  of  using  an  ideal  matched  filter  for  pulse  compression. 

The  right-most  peak  on  each  trace  corresponds  to  the  scattering  center 
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RANGE  IN  METERS 
(b)  TARQET  AT  2100  km 


FIG.  38  ALCOR  linear  video  ( J RCS  ) return  from  two  element  target. 


RANGE 


_(b)  TE 

FIG.  39  ALCOR  video  return  from  a cylinder  (average  S/N  80  dB 

_ ■*>  r"  . > 

-87-  • 


F/6  17/9 


most  distant  from  the  radar  while  the  peaks  on  the  left  are  the  responses 
from  the  two  scattering  centers  nearest  the  radar  on  the  end  of  the  cylinder. 

As  expected,  the  responses  from  these  scattering  centers  tend  to  coalesce 
near  axial  incidence  (top  trace).  At  this  viewing  angle,  the  range  profile 
extent  shows  the  length  of  the  target.  As  the  aspect  angle  is  increased 
to  approach  a broadside  view,  the  range  (or  line  of  sight)  extent  of  the 
target  becomes  indicative  of  the  diameter  of  the  cylinder. 

The  signal  to  noise  ratio  (SNR)  in  these  plots  is  determined  by 
the  polarization  and  aspect-dependent  radar  cross  section  (RCS)  of  each  of 
the  three  scattering  centers  as  well  as  the  distance  between  the  target 
and  ALCOR.  Thus,  to  specify  the  SNR  in  these  plots,  six  aspect  dependent 
numbers  are  needed.  Rather  than  quote  these  numbers  or  specify  the  SNR 
for  a 0 dBSM  sphere  at  the  same  range  as  was  done  in  the  previous  subsec- 
tion, we  provide  an  average  value  as  follows.  The  largest  RCS  shown  in 
Figure  39a  is  about  -11  dB  relative  to  one  square  meter  and  the  SNR  is 
69  dB.  A one  square  meter  sphere  at  the  same  range  would  have  been  an 
SNR  of  100  dB.  The  smallest  RCS  in  Figure  39a  is  about  -30  dBSM  and 
has  an  SNR  of  70  dB.  Thus,  the  average  SNR  in  Figure  39a  is  about  80  dB. 

The  rear  scattering  center  in  the  TE  case  (Figure  39b)  is  not  visible 
near  axial  incidence,  but  SNR  - 80  is  still  representative  of  the  level 
of  noise  in  this  picture. 

Figure  40  shows  a much  lower  SNR  (the  target  is  more  distant 
from  ALCOR) . The  largest  and  smallest  SNR's  in  Figure  40a  are  29  and 
10  dB  so  that  the  average  SNR  is  about  20  dB.  Figure  40b  shows  the  same 
situation  except  that  TE  polarization  is  used.  Note  that  the  rearmost 
scattering  center  is  barely  visible  at  this  distance  from  ALCOR. 
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SECTION  6 
SUMMARY 


In  this  study  the  technique  for  computing  the  smoothed  impulse 
response  has  been  extended  to  the  case  of  scattering  from  triangular 
plates,  quadrilateral  plates,  and  more  refined  aircraft  models  that  consist 
of  a cylindrical  fuselage  with  triangular  or  quadrilateral  plates  for  the 
wings,  stabilizers,  and  rudder.  To  the  knowledge  of  these  writers  this  is 
the  first  time  that  a successful  solution  of  these  problems  has  been  ob- 
tained. The  numerical  techniques  that  were  developed  consisted  of  the 
simultaneous  solution  of  two  vector  space-time  integral  equations.  A 
major  effort  was  required  in  order  to  be  able  to  handle  non-rectangular 
plate  geometries.  The  resulting  numerical  procedure  consisted  of  a set 
of  local  coordinate  transformations  in  the  plate  edge  regions  in  order  to 
be  able  to  properly  satisfy  the  boundary  conditions  on  the  surface  currents. 
Results  were  computed  for  several  triangular  plates,  finned  cylinders, 
and  a MIG-21  aircraft  model.  These  results  were  found  to  be  in  good 
agreement  with  measurements  taken  on  the  Sperry  time  domain  scattering 
range . 


At  high  frequencies  where  the  characteristic  target  dimensions 
are  many  wavelengths,  the  target  response  is  due  to  returns  from  only 
local  regions  of  the  target  surface.  These  regions  occur  where  the 
derivatives  of  the  surface  tangent  are  discontinuous  and  are  referred  to 
as  scattering  centers.  They,  in  turn,  produce  singularities  in  the  target 
impulse  response.  In  this  effort  time  domain  techniques  that  utilize 
the  local  nature,  both  in  space  and  time,  of  these  returns  were  developed 
and  demonstrated.  Both  the  numerical  solution  via  the  space-time  integral 
equation  and  the  measurement  solution  via  the  time  domain  scattering  range 
were  carried  out.  Results  were  obtained  for  several  test  target  scattering 
centers.  The  power  of  this  time  domain  approach  to  target  response 
prediction  can  not  be  overemphasized.  It  provides  an  efficient  analytical 
and  measurement  capability  for  obtaining  accurately  the  response  of  large 
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targets.  In  addition  measured  results  were  obtained  for  the  response  of 
scattering  centers  when  coated  with  RAM.  These  results  show  that  RAM 
can  indeed  lower  the  scattered  response  but  that  the  amount  of  reduction 
is  dependent  upon  the  spectrum  of  the  illumination  and  the  aspect  angle. 
Moreover,  there  were  several  frequencies  where  the  RCS  experienced  nulls 

I of  an  additional  8-10  dB. 

The  outstanding  achievement  of  this  work  has  been  the  development 
of  general  techniques  for  computing  the  wideband  (or  impulse)  response  of 
realistic  targets  with  a highly  efficient  and  accurate  method  of  computing 
the  noise-corrupted,  pulse-compressed  video  response  of  these  targets  when 
illuminated  by  a wideband  radar.  In  particular  the  application  of  some 
subtle  properties  of  the  discrete  Fourier  transform  enabled  the  calculation 
of  transforms  of  large  time-bandwidth  product  waveforms  (such  as  5000) 
using  an  FFT  with  less  than  one- tenth  the  number  of  points.  These  tech- 
niques were  demonstrated  using  a model  of  the  ALCOR  radar. 
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SECTION  8 
APPENDICES 


8.1  SELF- PATCH  INTEGRATION  AT  EDGE 


We  need  to  find  the  integral  over  the  self-patch 


t (r  t)  - -L  f fUMl 

As  re(t  Try  r 


, T « t - R 


(8-1) 


as  used  in  Section  2.3.2,  where  rg  is  a node  on  the  edge;  and  also 


-Jeff  .-*■ 
J (r 


= Ai "/ 

P ^ 


J(r ' ,T)  dS' 


(8-2) 


We  assume  a semicircular  patch  of  radius  p centered  at  r^  and  assign  a 

local  coordinate  system  £,  n,  II  and  i to  the  edge.  See  Figure  41. 

Note  that  As  = it/2  p2  or  p = / 2AS  /it  ' . There  is  a reference  node  at 
0 ^ V 6 

r , which  is  separated  from  r by  the  perpendicular  distance  d.  Also, 

° e -► 
dS'  - p'  dp'  d<(>  and  R = p' . In  general,  J cannot  be  assumed  to  be  con- 
stant over  the  patch,  but  must  satisfy  the  edge  conditions  (11) . Perform 
the  integration  by  components  ( II  or  i)  . 

(a)  J is  constant. 

This  condition  can  be  assumed  to  hold  for  at  a bound  edge, 
where  the  boundary  condition  is  3J^/3ri  * 0.  Then  simply: 


_eff  .-*•  . 

J (r  ,t) 
e 


(8-3) 


and 
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(8-4 


W*’  = 4?  J 


4 4%  IT  o 


where  J is  the  value  at  the  reference  point  (r  ,t) 
o ° 


(b)  J = 0 at  edge. 


This  condition  holds  for  at  a free  space  edge.  Assume 

J(n)  = h/d  J and  that  there  is  no  variation  of  J in  the  direction  t. 
o 

Using  r|  = p’  sin  4>: 


ef f -*■ 

J (r  ,t) 
e 


p * J / 3\  / \ 

f p’2  dp ' J sin  4>  d4>  = y2  J 


eff  -»■  . _ 4 

J (r  ,t)  = —7 
e 3ird 


Similarly 


VVt(  -35  / o'  *'  / sin  * d* 


(2 AS  /n) 

J 

4ird  o 


We  will  write  this  as 

, /2AS  1 „ _ 

Vv«  - Te  i^r 

[The  condition  J = 0 is  also  satisfied  by  the  ll  component  at  a bound  or 

ef  f 

symmetry  edge.  For  this  case  we  prefer  to  set  J =0  and  Ag  = 0,  since 
the  image  of  the  self  patch  has  currents  with  opposite  sign  in  the  parallel 
component.]  The  effective  current  is  close  to  zero  and  its  exact  value 
does  not  matter  much. 
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(c)  J ~ J—  at  edge. 


This  condition  holds  for  the  parallel  components  at  a free  edge 
Assuming  no  variation  in  the  £ direction 


Using  n = p'  sin  <j> 


ru/4) 

r(3/4) 


Similarly  obtain 


8.2  FINITE  ELEMENT  TECHNIQUE 


The  following  is  an  adaptation  of  the  methods  presented  in  [13] 


We  wish  to  solve  the  differential  equation 


on  the  surface  (here  contained  in  the  x,z  plane) , subject  to  boundary 
conditions  on  J,  where 


The  finite  element  method  proceeds  by  dividing  the  surface  into  elements 
of  arbitrary  size  and  shape  with  a certain  number  of  nodes  on  the  boundary 
of  the  element.  Here  we  will  assume  that  3 nodes  define  the  element,  i.e. 
the  elements  are  triangles.  It  is  readily  seen  that  each  node  on  the 
surface  usually  belongs  to  6 triangles  for  a more  or  less  regular  (e.g. 
hexagonal)  tessellation  of  the  surface.  The  use  of  3 nodes  per  element 
will  result  in  A continuous  to  first  derivatives  across  element  boundaries 
(For  higher  order  continuity,  more  nodes  per  element  would  need  to  be 
used. ) 


where  N ^ are  simple  interpolating  functions,  and  where  the  summation  is 
over  the  3 nodes  of  element  (e) . Note  that  time  and  space  variation  has 


been  separated  in  this  expression 


a.  + b.  x + c.  z 
.,e . , 1 i 1 

N.(x,z)  

i „ » e 


2A 


(8-12) 


a . = x . z,  - x z . 
i D k k ] 


b.  = z.  - z, 
ilk 


(8-13) 


C.  = X,  - X. 
i k ] 


2A®  = 


X.  Z . 
1 1 


X . z . 

3 D 


= twice  area  of  element 


(8-14) 


x.,  z. , etc.  are  the  coordinates  of  the  3 nodes. 

l l 


T0 

Substituting  (8-11)  into  (8-9)  yields  f(A  ) f 0 which  is  called  the  resi- 
dual. The  method  of  "weighted  residuals"  using  Galerkin  weighting  N^, 
then  assures  a stable  and  minimized-error  solution;  that  is,  we  are  to 


solve  for  A . such  that 
1 


/ 


N6  f(A®)  dS  = 0 
1 


(8-15) 


for  all  nodes  i of  element  (e) . We  consider  the  above  equations  first 

for  one  element,  then  the  equations  for  all  elements  will  be  combined 

-► 

to  give  one  large  matrix  expression  for  the  A^ . 

Noting  that  (8-9)  contains  only  surface  components,  substitute 

“►g 

f (A  ) in  (8-15).  For  the  x component,  we  obtain 
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e (V  • *6)  ds 


3 2 A* 
X 


f Ne  Y dS  + ( 

-/ e 1 3t  Je 


3E1 

ds  + / e Ne  -rr  ds  = D 


i 3t 


(8-16) 


We  now  apply  the  integration  by  parts  formula: 


y*  N*  nxVD  dS  = y*  (n  x n^)  ND  dl  - J"  D n x VN  dS  (8-17) 

S C S 

where  n (with  components  n and  n ) is  the  normal  to  the  contour  C of 
c x z 

region  S.  Equation  (8-16)  then  becomes 

3nS  / 3n6  3n®  \ 

/ N*  £ N.(V  • A*).  nx  d*  - y £ 17  Axj  + Azj/dS 
*'ce  j se  j 


3n 


3a 
’j  9y 


ym  wm  . 

1 

P 3x 


3Ne  3ng 


3n6  3n6 


/u«  . un  . 

s- 

:®  * f Ne  Ne  ds 

13  Je  1 3 

S 

ex  f »,e  „e 

).  . = I N.  N . n 

13  J*  1 3 x 


(8-19) 


y.  Bn1: 

— «T“  Ne 
3x  3 


Summing  over  all  elements  which  contain  both  i and  j , define 


k**  = V * 

13  !3 


(8-20) 


and  similarly  for  the  other  k's,  b's  and  t's. 


Hence,  summing  equations  (8-18)  for  all  elements  we  cam  write 


a matrix  equation 


MM  - MM  - Mkl  - * 4^1  ■ ° 


(8-21) 


where  k*^  is  an  element  of  etc.,  where  Jv  • a"|,  J^ax1'  etc.  are 

column  vectors  of  values  at  the  nodes.  There  is  a similar  equation  with 
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x and  z interchanged  for  the  other  component.  Note  that  in  performing 
the  summation  over  elements  the  contributions  from  the  contour  integrals 
cancel,  except  for  the  contributions  from  the  problem  boundary.  That  is, 
0X]  contains  only  entries  at  i,j  both  on  the  problem  boundary,  and 
V • A needs  to  be  evaluated  only  at  the  problem  boundary. 

Writing  the  time  derivative  as 


92A 


A (t  + At)  - 2A  (t)  + A (t-At) 
X . . _ _X X x 


(8-22) 


at 


(At)' 


combining  like  terms  and  multiplying  by  T \ obtain  the  matrix  equation: 
pMt  + Atrj  = -pWt-At)]  + £21  - (At)2  T_1  rXXJ  fAx  (t)"j  + (At)  2|e  (t)| 

- ^(At) 2 T_1  YX]f-^-(t)]  + [(At)2  T_1  BX]fv*A(t)] 


- [(At)2  T-1  KXZj[Az(t)| 


(8-23) 


and  a similar  equation  for  the  other  component. 

The  matrix  expressions  in  brackets  are  all  completely  defined 
by  the  geometry  and  hence  need  to  be  determined  only  once  at  the  beginning 
of  the  problem. 

This  is  the  solution  to  the  original  differential  equation. 

It  is  similar  to  the  finite  difference  solution  presented  in  Section  2, 
except  that  surface  differentiation  has  been  replaced  by  matrix  multipli- 
cation . In  other  respects  the  solution  is  the  same  as  in  Section  2 : the 
surface  current,  J,  at  all  interior  nodes  is  obtained  by 

‘ Y (*»S(t> 
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J(t) 


where 


\s 


(t) 


= _L  f lill^Uds' 

4ir  J R 


NS 


The  boundary  conditions  enter  into  (8-23)  via  the  values  of  (V  • A)  at  the 
boundary.  We  can  write  this  as 


V • A = V 


j_  r !ids,  =j,  f v . m ds. 

47T  J R 4TT  y \ R / 


(8-24) 


This  integral  is  separated  into  self  and  non-self  parts.  The  self  part 
is  evaluated  by  applying  the  boundary  conditions  on  J,  in  a manner  similar 
to  the  procedure  in  Section  2.  The  self-integral  is  difficult  when  the 
integrand  is  the  parallel  component  of  J at  a free  edge,  since  it  -*■  °°; 


it  can  be  evaluated  by  noting  that  + C,  9J||/9£  where  £ is  the 


coordinate  l|  to  the  edge,  and  by  noting  that  the  contribution  of  J 
zero. 


Ilo 


is 


This  solution,  just  as  the  one  in  Section  2 thus  requires 
special  treatment  at  the  edges.  Further,  it  requires  a large  matrix 
multiplication,  instead  of  simple  differencing,  at  each  time  step.  The 
advantage  is  that  any  shape  surface  can  be  treated  and  that  the  solution 
should  be  more  accurate  ( and  hopefully  more  stable)  due  to  the  fact  that 
A is  continuous  to  first  derivatives  over  the  surface. 

8.3  FIRST  ORDER  SCATTERING  CENTER  RESPONSE 

The  impulse  response  of  a target  in  the  physical  optics  approx- 
imation can  be  shown  to  be 


-102- 


(8-25) 


foh  (£*\  = -L  d2s 

• W ' 2TT  c2dt2 

where  S(t)  is  the  illuminated  area  projected  on  the  plane  normal  to  the 
direction  of  incidence.  The  above  expression  is  polarization  independent. 
By  considering  the  interaction  of  surface  currents  near  the  specular  point, 
one  obtains  a first  order  polarization-dependent  correction  to  the  above. 
The  first  order  impulse  response  of  a specular  point  on  a smooth  curved 
patch  is  derived  in  Ref.  1.  With  our  usual  normalizations  it  is  given  by 


a a 


, 


where  T is  the  time  of  incidence.  The  first  two  terms  are  the  physical 
z 

optics  terms,  the  last  is  the  first  order  correction?  the  - sign  refers 
to  TE  polarization,  the  + sign  to  TM  polarization  (see  Section  3.2  for 
coordinate  conventions) . 

For  a prolate  spheroid  with  a,b  the  semi-minor  and  -major  axes: 


S(t) 


it  a2  b 


( ct^B  ~ c2t2/4 ) 
„3/2 


u ( t) 


(8-27) 


where 

2 2 2 2 
B = b cos  a + a sin  a, 

u(t)  is  the  step  function, 

a = angle  of  incidence  with  respect  to  the  major  (z)  axis. 
CT  - ET 

— 5.  . Zil°.  ' is  the  time  of  incidence  for  the  target  centered 
a a at  the  origin. 
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Wt 


Evaluating  the  terms  in  (8-26) , we  obtain 


1 dS  I ab  . . _ . 

P = r — I ■ ^ (impulse) 

z 2TTa  cdt|  2B 
o 


1 d S 


,2b 


* 2’  c2dt2 


4B 


3/2 


(step) 


(8-28) 


.*  Sin2  . 


, 2,2 

2a  b 


z2 


^!o  ■ 4 “(-S-1) p-  ,s“p  cort*cti°"’ 


with 


f h (t)  ■ PZ  6141  * (*.  4 S«)  UtTI 


+ • • • 


For  reference  we  note  that  the  next  higher  order  term  in  the 
expression  (8-26)  is  a ramp.  The  first  order  correction  contributes  to 
this  term  (see  Ref.  1) ; we  can  also  expect  a contribution  from  the  second 
order  correction.  We  will  write  this  term  as 


(R  - + R 


z2  z3 


j T U 


(T)  , 


(8-29) 


,3  2 /K1  ~ *2 

Rz2  “ 2 Pz  3 \ 2 


i— - — -j  ± a|— — -js^  (ran?)  correction) 


R is  to  be  obtained  from  the  second  order  correction. 
z3 
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/ 


We  will  now  derive  the  area  function  of  a right  circular  cylinder 
of  radius  a,  at  oblique  incidence  (0  < a < v/2) . It  is  readily  seen  that 
S(t)  is  a piece  of  an  ellipse  given  by 


(i)  * ( 


£ cos  a\ 


(8-30) 


where 


x is  the  usual  coordinate , 

£ is  the  coordinate  in  the  yz  plane. 

At  time  t'  after  incidence  on  corner,  the  ellipse  is  truncated  by  a ver- 
tical line  (in  the  end  of  the  cylinder)  at 


^l^fc  * - cos  a 2 sin  a cos  a ' 


r = « major  axis  of  ellipse. 

o cos  a 


Writing 


we  have 


x - a /I 


- ■ / 


Cx(t) 


2p  dC  > u(t') 
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(8-31) 


Hence 


dS  _ _dS 
dt' 


«.  f1 


n 

1 

i t’  ) 

O 

4a  sin  ct) 

dC1  dt 


3/2 

(sin  a)  cos  a 


u(t') 


d2s 


dt 


,2 


rg  t'"^  (1 — — ) 

/ a Z V 8a  sin  a / 
3/2 

2 (sin  a)  cos  a 


u(t') 


(8-32) 


(8-33) 


We  thus  obtain  for  the  normalized  impulse  response 


— h (— ) 

= 1 

(» 

+ H 

a V a / 

r 

2 

cT 

ct 

z 

T 

a 

a 

u(T) 


(8-34) 


z 3/2 

4Ti(sin  a)  cos  a 


K = 
z 


5/2 

32TT(sin  a)  cos  a 


H , K are  to  be  obtained  from  the  first  order  correction  (not  performed 
z2  z2 

here),  and  contain  the  polarization  dependence.  For  a cylinder  of  length 
L and  radius  a centered  at  the  origin 


— — * - | — cos  oi+2  sin  aj 


(8-35) 


At  broadside  incidence  (O  = IT/2)  we  have 


S(t')  =2 


(8-36) 


-106- 


Differentiating  this  twice  we  obtain  the  polarization  independent  terns 


G - 

z 


L/ii 

2V 


H 

z 


3 L 

32v  /I 


where  the  response  is 


H , is  to  be  obtained  fron  the  first  order  correction. 
z2 


At  axial  incidence : 


S(t')  = x a2  u(t') 

2 

-d  = it  a2  6’(t')  (a  doublet) 
dt’ 

We  can  expect  the  first  order  correction  to  produce  the  next  higher 
order  tern  (impulse)  and  also  terns  of  the  fom  T and  higher  due  to 
the  edges.  That  is,  at  axial  incidence: 


— h ( — \ m D <S * ( T)  + P S(T)  + H , u(T)  + 


• • • 9 


